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correlated  with  particle  energy  input  to  within  £  3  sec,  and 
the  yield  varies  between  0.6%  and  1.3%.  No  statistically  sig¬ 
nificant  correlation  was  found  between  this  SWIR  emission 
efficiency  or  delay  and  a  series  of  auroral  (and  measurement) 
parameters,  which  included  the  fractional  solar  illumination  of 
the  radiating  air  volumes.  The  distributions  of  visible 
radiance  from  EXCEDE  are  found  to  differ  from  predictions  of 
accepted  models  of  electron  transport  and  atmospheric  excitation 
both  in  individual  photographic  frames  and  as  a  function  of 
charge-ejection  al t i tude;  pr inci pal ly  ,  there  is  an  excess  of 
radiation  from  within  a  few  rocket  diameters,  and  the  mean 
brightness  changes  much  more  slowly  than  the  ambient  air  density 
between  85  and  123  km.  Space-  and  t ime-dependences  of  the  01 
6300  A  glow  that  resulted  when  ionosphere- recombining  molecules 
were  exhausted  by  HEAO-C's  launch  rocket  (Eastern  Test  Range, 

20  Sep  79)  were  derived  from  calibrated  video  images;  its 
lateral  growth  rates  indicate  that  the  widespread  F-region  after 
glow  results  principally  from  reactions  of  H2>  Some  variability 
in  the  lower  thermosphere's  density  structure  is  inferred  from 
a  statistical  analysis  of  the  solar  photon-excited  vacuum  ultra¬ 
violet  radiation  intensities  measured  from  the  low  earth  orbiting 
satellite  S3-**.  Analysis  of  a  calibrated  ground  photograph  of 
the  visible  afterglows  from  PRECEDE's  pulsed  electron  beam 
unexpectedly  shows  the  energy  yield  to  be  at  least  comparable  to 
that  in  prompt  air  fluorescence  (^3%). 
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The  objective  of  the  work  reported  here  is  to  evaluate  and  inter¬ 
pret  data  taken  from  rockets  and  aircraft  on  the  IR-optical  radiations 
resulting  from  natural  and  artificial  energy  inputs  that  simulate  input 
from  atmospheric  nuclear  explosions.  Data  from  the  EXCEDE:  Spectral  and 
PRECEDE  rockets,  the  1 97 9- 1 980  aurora  and  a i rglow  mi ss i ons  of  AFGL's 
opt i cal ly- i nst rumented  aircraft,  and  USAF  satellite  S3_i*  are  addressed. 

The  field  measurements,  sponsored  by  the  Defense  Nuclear  Agency  and 
directed  by  Air  Force  Geophysics  Laboratory,  are  described  in  earlier 
DNA  reports  (Ref's  1 ,  21 ,  22,  27,  30,  and  Al). 

Measurements  of  auroral  particle-associated  2.8  -  3.1ym  enhancements 
with  6  mi  1  1  i  rad-field  radiometers  operated  at  aircraft  altitude  (M2  km) 
show  that  l)  the  SWIR  output  is  correlated  with  column  energy  input  to 
within  .$  3  sec  and  2)  the  mean  energy  efficiencies  vary  between  0.6  and 
1.3%  (a  similar  factor-2  variation,  in  a  lower  range,  was  observed  in 
earlier  aircraft-based  measurements).  No  statistically  significant 
correlation  was  found  between  SWIR  yield  or  "lag"  and  auroral  or  measure¬ 
ment  paiameters  such  as  fractional  solar  illumination  of  the  emitting 
region,  intensity  of  particle  precipitation,  type  of  auroral  form,  and 
instrument  pointing  direction  (along  or  at  angles  up  to  26°  from  the 
geomagnetic  field  lines).  The  aircraft  data  provide  some  indication  of 
an  increase  in  SWIR  emission  efficiency  with  energy  deposition  altitude. 

The  visible  radiance  distributions  in  air  excited  by  EXCEDE:  Spectral's 
3  keV  electron  beam  are  found  to  differ  from  predictions  of  an  independ¬ 
ent  particle  transport  model  both  within  individual  onboard-photographic 
frames  (excess  radiation  close  to  the  rocket)  and  as  a  function  of  injection 
altitude  (anomalously  narrow  range  of  mean  surface  brightnesses).  Spatial 
distributions  of  volume  emission  rate,  unfolded  from  the  isophote  plots, 
can  serve  in  determining  the  glow's  excitation  mechanisms  --  discharge- 
or  simple  outgass i ng-rel ated  --  and  thus  in  assessing  the  contribution 
of  secondary  electrons  to  the  measured  infrared  excitation. 
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Low  light  level  video  image  of  the  01  6300  A  afterglow  produced  by 
ionosphere  recombination  reactions  following  release  into  the  F  region 
of  rlj  and  H^O  vapor  by  the  20  Sep  79  HEAO-C  launch  rocket  were  converted 
to  radiance  contour  plots,  from  which  time-  and  space-dependences  of 
emission  intensity  were  derived.  The  lateral  growth  of  the  plasma- 
depletion  glow  suggests  that  it  is  caused  principally  by  the  more  rapidly- 
diffusing  molecules. 

An  analysis  of  the  Lyman-B i rge-Hopf i eld  band  intensities  from  the 
daytime  thermosphere,  measured  by  a  nadir-pointing  narrow-field  photo¬ 
meter  on  polar-orbiting  satellite  53-^,  shows  the  data  samples  to  have 
somewhat  higher  variance  than  would  be  expected  from  photon  counting 
statistics  alone.  This  finding  indicates  some  variability  over  distances 
25  km  in  density  of  the  emitting  layer  (135  "  1 80  km)  at  latitudes 
near  70°  N. 

h  groundbased  photograph  showing  visible  radiation  pers i st i ng  M 0 1  s 
sec  after  excitation  of  the  atmosphere  by  the  2^  keV  electron  pulses  from 
PRECEDE  was  quantified.  The  t ime- i ntegrated  light  yield  was  found  to  be  com¬ 
parable  to  or  greater  than  the  yield  in  prompt  fluorescent  radiations. 

The  source  of  this  long-lived  radiation  has  not  been  identified. 
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PREFACE 


This  is  a  report  of  the  first  year's  work  on  Contract  DNA001-81-C 
0003,  covering  evaluation  and  interpretation  of  nuclear  effects-simulav ing 
infrared-optical  data  taken  by  Air  Force  Geophysics  laboratory  under 
sponsorship  of  the  Defense  Nuclear  Agency.  The  information  comes  from  the 
following  field  experiments. 

—  EXCEDE:  Spectral  electron-inject  ion  rocket  EX8 5 1  - -'+4- 1 
(19  Oct  79) 

--  PRECEDE  electron- inject  ion  rocket  EX407 . ^ 1  - 1  (17  Oct  7A) 

--  USAF  NKC-135A  Aircraft  55“ 33 1 20  auroral  and  airglow 
missions  (in  1 979~80) . 

Atn. -spheric  radiance  data  from  experiment  CRL-726  on  USAF  satellite  S  3_/* 
(1978)  were  also  considered.  Preliminary  assessments  of  the  results  from 
EXCEDE  and  the  aircraft  are  reported  in  Refs  1  and  22. 

The  work  was  done  under  the  direction  of  i.L.  Kofsky.  Mrs.  C.C. 
Rice  typed  the  manuscript.  Critical  information  about  the  field  measure¬ 
ments  was  provided  by  E.R.  Huppi  and  R.R.  O'Neil  of  AFGL's  Optical  Physics 
(OPR)  Branch,  and  the  authors  benefited  from  discussions  with  D.N.  Anderson 
and  R.E.  Huffman  of  AFGL  and  M.  Mendillo  and  J.  Baumgardner  of  Boston 
Uni versi ty. The  encouragement  of  A.T.  Stair  Jr  of  AFGL,  and  Capt.  P.  Lunn , 
Lt.  Col.  W.  McKechney,  and  Dr.  H.C.  Fitz  Jr  of  the  Defense  Nuclear 
Agency  is  gratefully  acknowledged. 
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SECTION  1 


SPATIAL  DISTRIBUTION  OF  THE  HEAO-C  DEPLETION  AIRGLOW 

INTRODUCTION 

PhotoMetr ics ,  Inc.,  working  under  AFGL  direction,  made  a  series  of 
measurements  from  AFGL/DNA  NKC-135A  aircraft  5 5“ 3 1 20  on  the  airglows  that 
resulted  from  explosive  release  of  and  H^O  vapor  into  the  F  region 
by  the  Atlas-Centaur  rocket  that  launched  NASA's  High  Energy  Astrophys ica 1 
Observatory-C  satellite  in  the  early  morning  of  20  September  1979.  A 
preliminary  report  of  our  results  along  with  an  outline  of  the  ionosphere- 
diagnostics  program  is  presented  in  Ref  1,  and  data  from  other  parti¬ 
cipating  groups  are  reported  in  Ref  2.  The  spatial  distribution  of 
radiance  in  the  01  63 00  A  (forbidden  red)  line,  which  we  recorded  with 
a  low  light  level  all-sky  video  camera,  is  an  observable  of  particular 
usefulness  in  verifying  theoretical  models  of  the  plasma  density  de¬ 
pletions  that  are  produced  by  such  injections  of  reactive  molecules, 
several  more  of  which  are  to  be  performed  in  the  near  future.  In  this 
Section,  we  reduce  these  calibrated  images  to  "maps"  of  this  recombina¬ 
tion  radiation  and  derive  glow  growth  parameters  from  these  plots,  for 
the  purpose  of  validating  and  determining  input  parameters  for  computer 
calculations  of  ionosphere  modification  phenomenology. 

No  evidence  of  irregularities  in  electron  density  of  spatial  scales 
that  would  cause  phase  and  amplitude  scintillation,  in  UHF-VHF  satellite 
communication  signals  has  been  derived  from  the  optical  and  electro¬ 
magnetic  wave  propagation  data  on  the  HEAO-C  ionospheric  "hole,"  which 
is  a  factor-rv5_deplet ion  region  surrounding  the  launcher's  ^2000  km  long 
trajectory.  Theory  (Ref  3)  indicates  that  fluid  instabilities  will 
develop  when  the  energy  source  perpendicular  to  the  ion-confining  geo¬ 
magnetic  field  (neutral  wind  for  gradient-drift,  ionospheric  currents  for 
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current-convective,  and  gravity  for  g rav i tat iona 1  Rayleigh-Taylor  in¬ 
stability)  is  sufficient  to  overcome  the  quenching  effect  of  free- 
electron  recombination  by  the  polyatomic  ions  that  result  from  release 
of  reactive  molecules.  In  consequence  the  Air  Force  Geophysics  Laboratory 
will  shortly  attempt  to  produce  artificial  equatorial  spread-F  irregular¬ 
ity  with  two  explosive  releases  of  H^O  and  CO^  vapor  from  Natal  (the 
Brazil  Ionospheric  Modification  Experiments,  currently  planned  for 
September  1982).  Additionally,  during  the  Spacelab  2  mission  (Fall  1984) 
Shuttle  Orbiter  is  scheduled  to  inject  reactive  gases  over  five  iono¬ 
spheric-observation  sites,  by  firing  its  Orbital  Maneuvering  System 
engines.  Reaction-rate  and  atmosphere-dynamics  parameters  are  needed 
for  predictive  models  of  these  ionosphere  modifications. 

BACKGROUND 

For  orientation,  we  review  our  optical  data  and  the  other  experi¬ 
ment  information  on  the  transport  of  plasma-recombining  molecules  from 
HEAO-C's  launch  rocket. 

The  aircraft  flew  parallel  to  and  approximately  40  km  N  of  the 
"midpoint"  of  the  rocket  exhaust  track  in  the  first  hr  after  launch 
(see  Fig33of  Ref  1).  Its  monochromatic  video  camera  (Table  1  and  Ref  4) 
recorded  moderately  high  signal/noise  6300  A  images  from  min  after 
initiation  of  second-stage  fuel  burn  (at  212  km  altitude)  until  at  least 
one  hr  after  the  vapor  had  been  injected.  The  red-line  afterglow  ex¬ 
panded  to  fill  the  all-sky  field  (practically  speaking,  ''-•1000  km  hori¬ 
zontally  in  each  direction),  developing  a  sha 1 1 ow  m i n i mum  roughly 
centered  on  the  trajectory.  A  2°  circular  field  wavelength-cycling 
photometer  pointed  to  the  zenith  measured  maximum  enhancements  of  8.5 
ki lorayleighs  in  the  01  forbidden  red  line  90  ±  10  sec  after  the  rocket 
passed  ^35  km  from  its  field  of  view  and  0.95  kR  in  the  green  (5577  A)  line 
50  ±  12  sec  after;  both  radiance  pulses  had  duration  halfwidths  near  2  min. 
Fig  1  shows  the  radiances  in  the  aircraft's  zenith  at  the  onset  of  the 
recombination  glows;  the  data  for  the  .  ull  flight  are  in  Fig  37  of  Ref  1, 
at  a  more  compressed  time  scale. 
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Table  1.  All-sky  video  camera  specifications 


Type: 

Dynamic  Range: 

Spectral  Response: 

Sens i t i v i ty: 

Optics: 

Geometric  Distortion: 
Recording: 


Cohu  (San  Diego,  CA)  2856B  ISIT,  16  mm 
ta  rget  (RCA  A8A9  A) . 

Automatic  gain  control  in  camera,  automatic 
brightness  cgntrol  in  intensifier  combine 
to  factor  10  . 

S-20  photocathode  preceded  by  narrow  band 
interference  filter  at  telecentric  stop 
of  the  camera  lens.  Center  sensitivity 
wavelength  6300  A,  bandwidth  ^20  A. 

0.1  kR-sec  with  5  £p/mm  at  photocathode 
(see  Ref  **)  . 

F/1.5  (T1.8  at  6300  A)  Pacific  Optical 
fish-eye  objective,  3.^5  mm  f£,  F0V  165°. 
Relay  lens  as  in  Ref  A. 

Maximum  within  a  circle  whose  diameter 
does  not  exceed  the  picture  height. 

Standard  analog  video,  1  inch  tape. 
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Figure  1 


Zenith  radiance  at  the  aircraft  near  onset  of  the 
enhancement  of  the  two  01  forbidden  lines.  The 
rocket  passed  about  30  km  S  of  the  photometer's 
field  of  view  at  the  time  marked  p. 
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The  aircraft  photometer  and  near-infrared  radiometer  data  so  far 
show  convincing  evidence  of  sky  background  enhancement  in  only  one  other 
wavelength  band,  2.832  -  3.14pm  (FWHM).  Several  other  airglow  features 

monitored  —  0,  (infrared  atmospheric  fundamental),  NO  +  0  continuum,  OH 

^  2  A 

av  =  4  and  5  sequences,  Na  I  resonance  doublet,  N  D  -  S  F-region  doublet, 

and  in  particular  the  (l 0 ,4)  hydroxy i  band  near  6700  A  reported  (Ref  5) 
as  excited  by  the  "Lagopedo"  (1977)  release  of  water  vapor  into  the 
ionosphere  --  were  not  detectably  increased  (as  far  as  can  be  deter¬ 
mined  from  the  present  stage  of  data  reduction).  We  interpreted  the 
^20  kR  2.832-3. 14pm  increase  (Ref  l),  which  was  measured  by  a  cryocooled 
radiometer  with  0.4°  field  pointed  13°  forward  of  the  aircraft’s  zenith, 

as  Av  =  1  hydroxyl  vibrational  emission  following  electronic  excitation 

2  + 

and  radiative  relaxation  of  the  OH  A(  Z  )  state.  This  initial  electronic 
excitation  presumably  results  from  the  e  +  H20+  dissociative  recombina¬ 
tion  reaction.  As  explained  in  Ref  1,  it  would  provide  a  selective 
population  of  the  upper  states  of  the  i  nf  rared  (1  , 0)  and  ( 2 , 1)  v  i  brat  i  ona  1 
bands,  to  which  the  radiometer  is  sensitive.  Our  interpretation  of  the 
emission  process  can  be  verified  by  spectrophotometry  of  the  A-X  bands 
of  OH,  some  of  which  lie  at  wavelengths  just  above  the  atmospheric  ozone 
cutoff  near  3100  A,  following  injection  of  H^O  vapor  or  into  the 
F  region.  (This  infrared-backgrounds  issue  is  not  further  pursued 
in  this  Report.) 

The  most  extensive  measurements  of  the  HEA0-C  plasma  depletion 
were  made  with  a  network  of  pol a r imeters  ,  which  monitor  by  Faraday 
rotation  of  radiowave  polarization  the  changes  in  total  free  electron 
column-content  along  lines  of  sight  to  beacons  on  geosynchronous 
satel 1 i tes .  These  proj ect ions  of  the  volume  electron  densities  to 
spatially-dispersed  stations  compose  low  resolution,  essentially  one¬ 
dimensional  "images"  of  a  limited  section  of  the  ionosphere  depletion 
region.  Similarly,  polarimetry  data  taken  using  earth-orbiting  satellite 
beacons  provide  1-D  scans  of  total  electron  content  over  the  limited 
time  span  of  the  satellite's  traverse.  (Additionally,  some  plasma 
depletions  near  the  F-region  peak  were  measured  directly  by  a  distant 
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incoherent-scatter  radar.)  The  monochromatic  all-sky  camera  views  of 
the  airglow  from  below  the  deposition  track  represent,  on  the  other 
hand,  high  spatial  resolution,  2-0  projections  of  the  column  concentra¬ 
tions  of  O^D,  which  are  a  by-product  and  indicator  of  the  removal  of 
ionospheric  plasma  by  the  rocket  engine  exhaust. 

THEORETICAL  BACKGROUND  AND  DATA  BASE 

Numerical  models  for  simulating  the  propagation  of  high-density 

chemica 1 1 y-react i ve  gases  injected  into  the  ionosphere  are  presented  in 

Ref's  6,7  and  8,  and  9.  These  (other  than  that  of  Ref  9,  which  was 

found  to  fit  the  data  poorly)  have  not  yet  been  applied  to  HEAO-C's 

initial  conditions,  which  result  from  what  is  in  effect  a  continuous 

explosion  along  the  boost  vehicle's  trajectory.  In  later  stages  of  the 

fuel  burn  the  net  exhaust  velocity  of  the  thruster  gas  may  be  in  the 

rocket's  direction;  and  some  major  (but  unknown)  fraction  of  the  ejected 

23 

H^O  vapor  condenses  into  droplets.  Some  3  x  10  molecules  were  in¬ 
jected  per  meter  of  fl ight  path. 

This  amount  of  material  makes  the  foreign  species  densities  com¬ 
parable  to  natural  ambient  in  the  MO's  km  radius  volume  to  which  the 
fast-moving  gases  expand  before  they  col  1 i s iona 1 1 y  thermal ize.  (Refer 
to  Figure  2  for  atmosphere  densities  at  trajectory  altitudes.)  A 
transition  to  diffusive  expansion,  which  is  more  readily  described  by 
transport  models,  takes  place  when  their  number  density  falls  signi¬ 
ficantly  below  ambient.  Ref  10  gives  the  theory  of  this  latter  growth 
phase  for  the  simplified  conditions  of  negligible  depletion  of  the 
atmospheric  reactants  and  diffusion  lengths  less  than  the  neutral- 
particle  scale  height,  and  applies  it  to  radiometric  photographs  of 
TMA/A50  vertical  trail  releases;  this  reference  provides  an  introduction 
to  both  the  physics  of  chem ica 1 1 y- react i ng  ,  cyl indrical ly-symmetric 
contaminant  clouds  and  the  data  access  and  reduction  procedures  applied 
here. 

The  and  H^O  molecules  destroy  both  themselves  and  the  iono- 
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sphere's  (predominant)  0+  ions.  (H^  is  also  removed  in  a  neutral- 
particle  reaction,  with  0  atoms  (Ref  11).)  At  the  same  time  the  0+ 
concentrations  in  the  "hole"  are  partially  replenished  by  diffusion, 
and  the  plasma  density  distribution  in  the  affected  volume  is  modified 
by  F-region  ion  and  neutral  winds  (Ref's  8,12). 

Fig  3&of  Ref  1  presents  sample  all-sky  images  of  the  6300  A-line 
glow  that  results  from  electron  recombination  of  oxygen-conta i n i ng 
molecular  ions.  The  geometry  of  these  video  photographs  is  illustrated 
in  the  side  and  plan  views  of  the  rocket  and  aircraft  trajectories. 

Figures  2  and  3*  The  images  are  as  usual  two-d i mens iona 1  projections 
to  the  camera  lens  of  volume  emission  rates  in  columns  extending  through 
the  three-dimensional  radiating  cloud.  (The  glow  volume  is  definitely 
optically  thin  to  its  own  radiation.)  Scene  radiances  tend  to  increase 
with  zenith  angle  from  the  aircraft  because  of  the  increasing  sight  path 
through  the  glow. 

Simple  van  Rhijn  corrections  of  the  type  used  on  natural  air- 
glow  (Ref  13)  can  be  applied  to  assess  the  vert ica 1 -col umn  radiances, 
such  as  would  be  measured  from,  say,  a  distant  satellite  in  the  glow's 
zenith.  The  release  trajectory  is  tilted  upward  rather  than  being  con¬ 
centric  with  the  earth's  surface,  and  is  such  that  the  the  amount  of  material 
being  exhausted  per  horizontal  km  is  decreasing  along  the  trajectory;  the 
cosines  of  the  track's  inclination  angle  are  shown  in  Fig  2.  Zenith 
(or  nadir)  radiances  can  be  found  directly  by  multiplying  the  projected 
radiances  by  the  cosine  of  the  zenith  angle,  which  is  the  angle  between 
the  line  to  the  camera  and  the  local  vertical;  when  radiances  from 
release  segments  at  different  altitudes  are  i ntercompared ,  the  afore¬ 
mentioned  correction  for  the  nonuniform  horizontal  exhaust  rate  must  be 
applied.  The  glow  radiances  can  thus  be  straightforwardly  converted  to  plan 
view  values  from  the  equi -brightness  plots  presented  here,  or  alterna¬ 
tively,  volume  emission  rates  calculated  from  the  ionosphere  modifica¬ 
tion  models  can  be  projected  to  the  camera  to  simulate  the  actual 
measurement  conditions. 
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(The  camera  was  diverted  for  MO  -sec  periods  to  image  the  glow  in 
5577  A  radiation  at  early  times,  and  to  image  at  6300  A  with  a  23°  x  ?7° 
field  at  later  t  imes .  These  br ief  data  segments  are  not  considered  here.) 

Horizontal-transport  information  has  been  extracted  from  the  video 
images  for  the  four  representative  exhaust  altitudes  257,  334,  400,  and 
450  km.  These  positions  along  the  launcher's  track  are  labeled  I  - 1 V 
on  Figures  2  and  3-  We  have  not  addressed  lower  trajectory  altitudes 
because  the  photogrammetric  accuracy  of  the  all-sky  camera  is  decreasing 
rapidly  near  the  edge  of  the  camera's  field  of  view;  as  Fig  2  shows,  the 
start  of  rocket-burn  at  212  km  altitude  lies  somewhat  outside  the  ^80° 
circular  field.  (Furthermore  the  6300  A  radiation  intensity  at  the 
lower  altitudes  tends  to  be  weak  because  0^D  is  both  strongly  collision- 
ally  quenched  and  unfavored  in  the  depletion  reactions,  as  we  show  in 
the  following  subsection.)  Near  exhaust  cutoff  the  transport  geometry 
loses  its  "cylindrical  symmetry,"  complicating  model  calculations.  The 
highest  point  we  considered,  450  km, is  950  km  hor i zontal 1 y  from  but onl y  8/1 0 
of  an  atmospheric  density  scale  height  below  the  end  of  the  rocket  fuel 
burn  at  501  km  altitude.  The  actual  mean  altitudes  to  which  the  trans¬ 
port  results  refer  are  higher  than  those  of  the  exhaust,  as  there  is 
net  upward  transport  of  reactive  gas  toward  lower  ambient  densities. 

The  four  ejection  a  1 1 i tudes span  an  atmospheric  density  range  of 
a  factor  31  (MSIS  II  atmospheric  model.  Ref  14)  and  an  ionospheric  [0+] 
range  of  a  factor  25  (Millstone  Hill  incoherent-scatter  radar  data, 

Ref  2).  ,  total  particle  (principally,  0  atoms),  and  0+  densities  are 

shown  at  each  data-reduct ion  location,  and  scale  heights  of  Nj  and 

H2°  (which  figure  in  the  diffusion  theory)  are  also  indicated  on  Fig  2. 

We  calculated  the  atmospheric  profile  at  the  time  and  latitude  of  rocket 
launch  using  MSIS  II;  results  are  listed  in  Table  2.  The  N^  molecule 
concentrations  were  used  in  calculating  the  probability  of  6300  A  radia¬ 
tive  relaxation  of  the  0  0  atoms  at  each  ^^20  exhaust  altitude;  since 
the  centroid  of  the  glow  is  above  these  altitudes,  there  is  actually 
somewhat  less  collisional  quenching  than  the^j's  indicate. 
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Also  shown  in  Fig  2  are  the  radial  diffusion  coefficients  of  H^O 
molecules  in  an  oxygen-atom  atmosphere,  as  given  in  Ref  12  for  1000  K 
exospheric  temperature  (we  derived  1025  K  from  MS  I S ) .  The  calculated 
diffusion  coefficients  of  are  closely  6  times  those  listed  for  H^O. 
Although  classical  minority  diffusion  through  a  uniform  atmosphere  does  not 
properly  describe  the  transport  of  the  injected  gas  (see,  for  example, 

Ref  15),  and  both  growth  while  it  is  initially  overdense  (Ref  1 7) and  deplet  ion 
of  the  critical  reactant  0+  must  be  considered,  the  diffusion  coefficients 
at  least  provide  some  orientation  to  the  order  of  magnitude  of  the  pro¬ 
jected  glow  dimensions.  For  example  at  10  min  after  exhaust  a  charac¬ 
teristic  spread  of  hydrogen  at  331*  km  (k  x  diffusion  coefficient  x  time)2 
would  be  about  200  km,  which  is  in  fact  close  to  that  calculated  for 
thermal ized  (see  Fig  7  of  Ref  8).  (The  cloud's  halfwidth  in  the 
minority-diffusion  phase  can  be  estimated  by  adding  in  quadrature  its 
halfwidth  at  the  end  of  the  self-diffusion  phase  and  the  halfwidth 
during  the  duration  of  minority  diffusion.) 

EXCITATION  OF  AND  RADIATION  FROM  0 1 0  ATOMS 

As  some  fraction  of  the  final-step  reactions  in  the  chains  that 
remove  the  dominant  F-region  ion  species  0+  results  in  production  of 
'o-state  oxygen  atoms,  the  rates  of  excitation  are  a  measure  of  the 
rates  of  ionospheric  plasma  depletion  (refer  to  the  discussion  in  Ref  1). 
Schematically,  considering  only  reactions  of  the  injected  molecules  and 
intermediate  species  with  ambient  charged  species 
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Considering  also  reactions  in  that  volume  close  to  the  trajectory  in 
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genation  and  electron  dissociative  recombination  of  molecular  ions 
(Kh/Kd,^10"2,  Ref  7), 


-2  1 

Under  these  conditions  --  10  [H_]  >,  [e]  --  excitation  of  0  D  is  not 

23  1 

favored.  We  note  that  10  H  molecules  are  exhausted  per  cm  of  trajectory; 

‘  -2  113 

when  these  expand  out  to  50  km,  10  [H.]~10  / nr  .  Comparing  this 

2  + 

figure  to  the  initial  ambient  electron  densities  (=  [0  ])  in  Fig  2,  we 
find  that  radiationless  plasma  depletion  becomes  important  within  a  few 
10's  km  of  the  release  trajectory,  particularly  at  the  lower  altitudes 
and  after  the  electron  density  has  been  reduced.  The  aforementioned 
fraction  (0  D  atom  excited  per  ion  removed)  thus  increases  away  from 
the  trajectory,  and  depends  on  the  F-region's  electron  density  profile. 

The  branching  ratios  for  excitation  of  0's  two  metastable  ground- 
configuration  states  in  the  electron  dissociative  recombination  steps 
are  not  known.  Pre-HEAO-C  theoretical  .  rguments  based  on  photodissocia¬ 
tion  data  and  channel  statistical  weights  (Ref  16)  indicated  that  M0- 
15^  of  these  reactions  excite  o'o;  however  the  actual  radiances  in  the 
red  line  were  considerably  higher  than  predicted  (in  Ref  9)-  Note  that 
the  ground  state  of  the  transition  that  produces  the  01  green  line 
emission  is  the  upper  state  of  the  red  line;  but  since  the  enhancement 
intensity  of  the  green  line  from  HEA0-C  was  a  factor  10  less  than  that 
of  the  red  line  (in  the  aircraft's  zenith.  Fig  1),  this  cascade  contri¬ 
bution  is  expected  to  be  small.  The  partial  rate  coefficients,  and  the 
impact  of  "close-in"  hydration,  can  be  estimated  by  fitting  parameters 
in  the  model  of  the  ionosphere  depletion  to  achieve  agreement  with  the 
o'd  surface  radiance  distributions  presented  here. 


The  factor-six  difference  in  diffusion  coefficients  of  the  two 
injected  species  might  suggest  that  radiance  traces  transverse  to  the 
launcher  trajectory  would  show  a  narrow  peak  superposed  on  a  broader  base, 
which  would  provide  a  more  or  less  direct  measure  of  the  relative  0 1 D 
excitation  probabilities  from  the  two  species'  reaction  chains.  On  the 
other  hand  the  explosive  nature  of  the  injection,  the  complicating  effects 
of  the  higher-than-ambient  exhaust  gas  densities  (which  include  the  afore¬ 
mentioned  successive  hydrogen  reactions  leading  to  Hi)"1"),  and  the  rapid 
reduction  of  [0  ]  among  other  factors)  may  wash  out  this  structure; 
indeed,  no  such  distribution  is  found  in  a  qualitative  examination  of 
the  video  images. 

The  metastable  oxygen  atoms  produced  in  the  mol ecul ar- i on  dissocia- 

1  3 

tion  reactions  radiate  at  6300. 3  A  (  0,  -  P  transition)  or  6363.8  A 

13  ^  ^ 

(  02  -  P-f)  >  or  undergo  collisional  quenching  by  N2  molecules.  (Deactiva¬ 
tion  by  02  can  be  neglected.)  Calculated  transition  probabilities  for 
the  two  red-line  radiations  are  O.OO69  and  0.0022  sec  1  respectively 

(Ref  18  )  ,  and  the  "aeronomic"  rate  coefficient  for  quenching  of  0  D  by 
- 1 1  2 

N2  is  3  ±  1  x  10  cm  /sec  (Ref  19;  the  uncertainty  is  probably 
larger  than  that  quoted,  as  the  several  measurements  vary  over  a  range 
of  a  factor  >  k) .  The  probability  that  O^D  atoms  radiate  in  the  video 
camera's  sensitivity  band  is  then 

_ 0.0069 _ 

0.0069  +  0.0022  +  3  x  10*11  [N2] 

where  [N^]  is  the  ambient  concentration  of  N2  molecules. 

The  calculated  6300  A  radiation  probabilities  at  injection  alti¬ 
tudes  are  listed  in  Figure  2.  Note  that  at  the  lowest  altitude,  212  km, 
only  1%  of  the  excited  atoms  radiate  in  the  camera's  sensitivity  band; 
at  257  km,  our  lowest  data  analysis  altitude,  i  radiate;  and  above 
^350  km,  which  includes  the  final  two  altitudes  considered,  the  fraction 
is  experimentally  indistinguishable  from  the  unquenched  figure  69/ 

(69  +  22)  =  0.76.  Qualitatively  speaking,  collisional  quenching  of  0  D 
reduces  the  glow's  radiance  only  at  zenith  angles  >  60°  in  the  direction 
of  the  start  of  the  burn,  and  (since  the  transport  of  injected  molecules 
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is  larqely  upward  and  horizontal  rather  than  downward  toward  increasing 
ambient  densities)  has  a  small  effect  on  intensities  over  most  of  the 
i mage  field. 

The  col umn-concentrat ion  of  0  D  atoms  is  given  by  (4"/0.0069) 

times  the  6300  A-line  radiance  in  the  direction  of  the  column  expressed 
2 

in  photons/sec  cm  sterad.  The  radiances  directly  indicate  production 
rate  of  the  end-product  species  at  the  (late)  times  that  the  fractional 
rates  of  change  of  [o'o]  due  to  both  exci tat ion-and-decay  and  mechanical 
transport  are  small  compared  to  1  in  the  characteristic  radiation  time, 
100  sec.  (The  lifetimes  of  OH  and  H^O  turn  out  ot  be  of  the  same  order 
when  the  plasma  concentration  [e]  is  near  its  ambient  value  of  lO'Vcm^.) 
Change  of  concentration  of  o'd  can  be  expressed  es 


d  [0  1  D  ] 
dt 


-  (0.0091  +  3  x  10-11  (N21)  to1 Dl 
+  (k][0H+]  +  [H  )  t  e ] 

+  .  •  v  [o’d]  , 


where  k^  and  k2  are  the  partial  rate  coefficients  for  excitation  of  the 

stave  in  dissociative  recombinations  (including  the  effect  of  cascading 

from  o's)  and  the  last  term  takes  into  account  the  net  transport  into 

Or  out  of  the  volume  considered  due  to  divergence  of  o'd's  flow  velocity 

v  (see,  for  example,  Ref  7).  The  0-atoms1  characteristic  diffusion 

7  2-1 

length  in  100  sec  is  of  the  order  of  50  km  (2  x  10  m  sec  x  100  sec)‘. 
Therefore  the  fractional  change  in  O^D  concentration  must  be  -  1  over 

this  horizontal  distance  if  the  excitation  rate  is  to  be  measured  directly 
by  the  radiative  decay  rate  (that  is,  when  excitation  and  decay  are 
proceeding  at  the  same  column  rate).  In  any  case,  the  metastable  state's 
long  lifetime  indicates  that  the  ionosphere  depletion  models  must  take 
into  account  transport  of  0^0  as  well  as  of  H^,  H^O,  and  intermediate 
ion  species  in  fitting  to  the  6300  A  radiance  data  distributions. 
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DATA  ACCESS  PROCEDURE 

Contour  plots  of  the  afterglow's  radiance  viewed  from  the  aircraft 
serve  as  the  starting  point  for  analysis  of  the  63OO  A-line  emission 
data.  Figure  b  i s  a  set  of  nine  equ i -rad i ance  plots  spanning  the  time 
period  from  the  first  evidence  of  a  glow  in  the  video  images  (125  sec 
after  second-stage  ignition,  385  sec  after  launch,  rocket  then  at  339  km) 
to  1060  sec  after  ignition  (780  sec  after  the  last  trajectory  altitude 
from  which  data  were  further  reduced  ).  Tabl e  3  lists  the  times  of  these 
plots  after  exhaust  injection  at  each  of  the  four  altitudes  considered. 

Note  that  the  aircraft  reverses  flight  direction  only  before  the  last 
point  in  the  data  set;  its  otherwise  straight  track  reduces  the  effort 
to  corre'  „  for  the  changing  view  perspective.  (Between  660  and  750  sec 
it  executed  a  25°  roll  to  starboard,  to  scan  the  zenith-pointing  radio¬ 
meter  transverse  to  the  rocket’s  trajectory;  only  one  image  frcm  this  j 

period  is  considered.) 

We  initially  concentrate  on  these  early  times  after  vapor  injection 
because  of  the  high  signal/noise  in  the  images  and  the  fact  that  virtually 

all  of  the  glow  remains  within  the  all-sky  (160°  circular)  camera  field.  j 

The  later-time  frames  also  have  useful  information  about  the  ionosphere  ! 

depletion,  particularly  because  the  production  and  decay  of  0 1 D  are  in  1 

near-equilibrium  and  the  aerochemi st ry  is  not  complicated  by  hydrogenation 

reactions  ([H^]  having  been  reduced  by  diffusion  and  reactions).  1 

The  brightness  maps  have  been  rad i omet r i ca 1 1 y  corrected  as  described 
here.  "Dimensions"  of  the  soft-edged  glow  and  effective  cross-track 
diffusion  coefficients  can  be  derived  from  cuts  across  the  2-D  plots, 
with  the  radiances  corrected  to  a  zenith  projection  to  remove  the  effect 
of  variable  pathlength  in  the  images.  (Attention  is  directed  to  the 

aforementioned  caveats  about  interpreting  the  afterglow’s  wid1’  terms  j 

of  species  transport  and  ionosphere  depletion.)  Further  descriptors  I 

of  the  space  and  time  dependence  of  the  recombination  qlow  can  of  course  I 

I 

be  taken  from  the  brightness  contour  plots  in  Fig  b. 
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Figure  4a -c  .  Contour  plots  of  the  6300  fl-line  afterglow  radiance  resulting 
from  HEAO-C's  launcher's  exhaust  of  H  and  H„0  into  the  iono¬ 
sphere,  projected  to  positions  1-3  of2 USAF  A7C *5S-3120.  4a 
projects  to  1  on  Fig  2 ,  4b  to  2 ,  4c  to  3.  The  radiances  are 
corrected  for  non-uniform  response  over  the  field  of  view 
(refer  to  text).  Absolute  radiancesare  the  product  of  the 
number  on  the  contour  lines  and  the  calibration  factor  given 
at  upper  right  of  each  plot.  The  nightglow  background  of 
130  -  350  R  (depending  on  zenith  angle)  has  not  been  subtracted. 
X's  indicate  the  straight  (horizontal)  lines  perpendicular  to 
the  launcher's  trajectory.  Refer  to  Fig  4d  for  a  zenith  angle 
scale. 
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Figure  **e. 


Contours  of  the  6300  A-line  afterglow 
projected  to  position  5  of  A/C  55“3120 
multiplier  is  2645R* 


radiance 
.  The 


27 


i 

i 

i 


Figure  4i.  Contour  plot  of  the  6300  A-line  afterglow  radiance 
projected  to  position  9  of  A/C  55*3120.  The 
multipl ier  is  AAO  R. 
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Table  3.  Times  of  (mage  radiance  contour  plots  (sec) 


Image 

T  after  L  T 

after  1 

1 1 

1  1 1 

IV 

1 

(339  -- 

o 

27 

50.5', 

67°  17.7'*) 

385 

85 

5 

- 

- 

2 

(400  -- 

27° 

51.1', 

67°  19.8') 

463 

163 

83 

3 

- 

3 

(437  -- 

27° 

53.7', 

67°  29.6') 

488 

188 

108 

28 

- 

4 

(«»84  -- 

27° 

55.6', 

67°  36.4') 

548 

248 

168 

88 

8 

5 

(480  -- 

27° 

57.2' , 

67°  42.4') 

6C0 

300 

220 

140 

60 

6 

(501  -- 

28° 

00.3'  , 

67°  50.0') 

700 

400 

320 

240 

160 

7 

( 

28° 

01.9'  , 

67°  59.7') 

750 

450 

370 

290 

210 

8 

( 

28° 

09.6'  , 

68°  28.6') 

1005 

705 

625 

545 

465 

[Ai 

rcraft  Turns 

180°] 

9 

( 

28° 

19.6'  , 

68°  05-5') 

1320 

1020 

940 

860 

780 

Altitude  of  At las-Centaur  rocket  in  km  *-  latitude  N,  longitude  W  of 
the  aircraft  measurement  platform. 


The  views  of  the  exhaust  trajectory  in  Figs.  2  and  3  indicate  these 
data-reduct ion  intercepts.  Horizontal  straight  lines  perpendicular  to 
the  track  are  curved  in  the  all-sky  projection.  Zenith  angles  for  Fig  4a-i 
a  re  s  hown  i n  F i g .  4a . 

This  first  reduction  of  the  video  camera  data  was  done  by  photo¬ 
graphic  photometry  from  a  display  monitor,  principally  because  thi s  procedure 
permitted  more  smoothing  of  the  "grainy"  scene  brightnesses  by  summing 
frames.  The  initial  analog  record  of  the  afterglow  has  also  been  copied 
onto  digital  tape  for  processing  using  AFGL/OPR's  video  image  analysis 
system,  whose  integration-storage  capacity  is  at  present  only  4  full 
frames  at  5“bit  radiometric  precision;  unless  image  data  are  stored 
offline  through  AFGL's  central  computer  (which  requires  further  pro¬ 
gramming),  this  analyzer  does  not  provide  sufficient  smoothing  over  a 
useful  dynamic  range.  In  practice,  30  full  video  frames  were  summed  by 
photographing  the  screen  for  1.0  sec. 
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The  data  tape  was  played  back  on  a  7s"  x  10"  (Tektronix  Model  632)  j 

i 

monitor,  with  display  gain  held  in  the  low  to  mid-range  to  maintain 
linearity  of  the  screen's  radiance  with  input  signal  voltage  (Ref  20). 

The  screen  was  photographed  from  2.3  m  with  a  long  focus  (200-mm)  lens  j 

so  as  to  minimize  radial  distortion  in  the  new  images.  A  fine-grain, 

high-contrast  35  mm  film  was  used  (EK  Pan-X).  Following  the  usual  calibra¬ 
tion  procedure,  we  placed  a  step  wedge  on  the  film  (also  at  1.0  sec  I 

exposure,  although  in  fact  this  precaution  against  error  due  to  recipro¬ 

city  failure  turns  out  not  to  be  necessary),  which  we  hand-developed  for 
uniformity  in  a  small  reel  tank.  The  density  distribution  on  the  35  mm 
images  were  accessed  with  our  two-dimensional  scanning  microdens i tometer , 
and  the  resulting  equidensity  plots  were  then  converted  to  distributions 

of  exposure  in  the  film-plane  using  this  step  wedge  calibration  (the  pro-  j 

cedure  is  the  same  as  that  for  the  EXCEDE:  Spectral  (Section  A)  and  earlier  |  J 

rad iomet r ica I  1 y-ca 1 i brated  photographs).  j 

Film-plane  irradiances  are  proportional  to  radiances  at  the  video  ’ 

| 

monitor,  which  to  the  extent  that  the  taping  and  playback  systems  respond  j 

linearly  are  proportional  to  the  initial  irradiances  at  the  video  camera's 

photocathode.  These  irradiances,  when  corrected  for  non-uniformity  of  | 

response  across  the  field  of  the  wide-angle  lens  and  taking  video  camera  J 

(and  also  the  display  monitor),  are  in  turn  proportional  to  the  radiances  i 

in  the  glow  cloud  --  the  phys ical 1 y- important  quantity.  We  made  this 
correction  (which  is  expected  to  be  same  in  all  frames,  that  is,  time¬ 
stationary)  by  adjusting  the  measured  natura 1 -background  6300  A  sky 
radiance  distribution  to  its  (van  Rhijn,  Ref  13)  dependence  on  zenith 
angle.  The  specific  v i gnet t i ng-cor rec t ion  procedure  is  as  follows.  In 
an  early  video  frame  (at  0535:48  UT ,  from  which  Fig  4b  was  derived, 
the  artificial  airglow  is  limited  to  a  small  (^10°  x  15°)  region  at  the 
(  western  edge  of  the  image;  the  remaining  exposure  is  due  to  the  natural 


6300  A  nightglow.  The  measured  apparent  irradiances  Z'  are  related  to 
the  irradiance  C'  at  the  frame  center  by 


where  V.  is  the  system  vignetting  or  overall  i nst rumenta 1  response (V.  ,<  1) 
and  V is  the  van  Rhijn  path-length  factor  through  the  nightglow  (V^  £  1), 

and  we  have  assumed  that  the  airglow  intensity  changes  little  over  the 
field  of  view  other  than  by  pathlength.  The  E^'s  were  read  point  by 
point  from  the  exposure  contours  in  the  test  frame,  and  the  van  Rhijn 
factors  were  taken  from  Ref  13  for  a  350  km  6300  A  emission  altitude 
(VD  varies  very  little  between  300  and  400  km,  the  expected  limits  of  the 
airglow  profile's  peak).  The  exposures  at  each  elevation-azimuth  posi¬ 
tion  in  the  data  frames  were  then  manually  corrected  by  the  factor 
E VVD(=  V.  C') .  The  maximum  value  of  V.  reaches  about  3,  near  some  80° 

K  i  i 

zenith  angles.  Note  that  this  correction  takes  into  account  the  two- 
dimensional  nonuniformity  of  response  by  the  imaging  and  reproduction 
systems . 

This  procedure  gives  relative  radiances  within  individual  video 
frames,  with  an  estimated  mean  accuracy  of  +  60%,  -  40%.  As  the  camera's 
automatic  gain  controls  were  used  to  allow  it  to  respond  to  the  necessarily 
wide  dynamic  range  of  scene  brightnesses,  an  absolute  calibration  for  each 
frame  was  taken  from  the  zenith  radiances  measured  by  the  2°-field  6300  A 
photometer  (shown  in  Fig  1  and  Fig  37of  Ref  1  ).  The  contours  in  Fig  k 
are  presented  as  successive  10  percent  decreases  in  scene  radiance  from 
the  radiance  at  the  aircraft's  zenith  when  the  maximum  radiance  was  near 
the  center  of  the  frame  (i.e.,  when  the  glow  comes  overhead),  and  10  per¬ 
cent  decrements  from  the  glow's  peak  in  earlier  frames. 

LATERAL  GROWTH  AND  TIME  DEPENDENCE  OF  GLOW  RADIANCE  AT  FIXED  POSITIONS 

The  i sophote  plots  show  the  expected  decrease  in  brightness  near 
the  launcher's  trajectory  at  the  lower  exhaust  altitudes,  presumably  due 
largely  to  depletion  of  the  ionosphere's  0+  ions.  They  also  indicate  a 
northward  drift  of  the  "axis"  of  the  afterglow,  which  we  had  initially 
ascribed  to  F-region  wind  (Ref  1).  Inconsistencies  in  the  direction  of 
the  "centerline"  of  the  glow  region  within  individual  frames,  however, 
lead  us  to  the  conclusion  that  this  apparent  motion  could  be  caused  by 


34 


errors  in  the  heading  and  other  attitude  angles  of  the  aircraft.  Its 
positions  were  taken  from  the  navigator's  log,  and  attitude  cannot  be 
checked  from  star  fields  as  few  stars  are  imaged  through  the  narrow 
wavelength-band  filter. 

Positions  along  the  one-dimensional  traces  (Fig  5)  have  been 
corrected  for  the  aircraft's  motion;  the  isophote  plots  in  Fig  4,  from 
which  the  traces  were  compiled, are  projections  to  the  individual  aircraft 
positions.  The  correction  of  elevation  and  azimuth  angles  was  made  with 
an  adaptation  of  a  navigation  program  designed  to  be  run  on  a  small 
digital  computer.  The  radiances  in  Fig  5  are  also  corrected  to  the  local 
zenith,  that  is,  the  variation  due  to  the  change  of  sight  path  length 
with  zenith  angle  has  been  removed  (as  discussed  earlier).  The  horizontal 
distance  scale  on  the  abscissa  was  assigned  with  the  assumption  that  the 
"centroid"  of  the  glow  in  the  vertical  direction  remains  at  the  initial 
exhaust  altitude  (as  also  noted,  a  net  upward  transport  is  expected). 

The  minimum  at  the  lower  exhaust  altitudes,  and  the  apparently- 
related  flattening  of  the  traces  at  higher  altitudes,  effectively  pre¬ 
cludes  direct  assignment  of  meaningful  halfwidths  to  the  radiance  traces. 
(Because  of  the  low  accuracy  at  257  km  exhaust  altitude,  which  is  due  to 
both  compression  of  the  zenith  angle  scale  and  uncertainty  in  the  vignetting 
correction  r  r  the  edge  of  the  all -sky  frame,  we  have  not  attempted  to 
derive  transport  data  from  these  lowest -al t i tude  radiances.)  We  used 
instead  the  separation  between  the  points  of  maximum  slope  in  Fig  5  to 
compute  effective  diffusion  coefficients.  (In  a  Gaussian  with  halfwidth 
to  1/e  maximum  /2a  --  i.e.,  described  by  exp-(x  /2a  )  --  the  slope  is  a 
maximum  at  x  =  a.)  The  slope  of  the  plot  o  tne  square  of  half  this 
distance  (i.e.,  the  half-width)  against  the  time  after  the  exhaust  mole¬ 
cules  are  deposited  can  be  interpreted  as  b  x  the  effective  diffusion 

coefficient,  as  described  in  Ref  10.  The  figure  derived  for  each  of  the 

7  2 

three  highest  exhaust  altitudes  is  about  2j  x  10  m  /sec,  which  is  close 
to  the  expected  diffusion  coefficient  of  near  400  km  (see  Fig  2). 

This  finding  suggests  strongly  that  the  wide  spread  of  the  glow  is  due 
principally  to  reactions  of  the  ejected  molecules  rather  than  to  the  H^O. 
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Figure  5a.  6300  A  local-zenith  brightness  along  lines 

perpendicular  to  the  direction  of  the 
launcher's  trajectory  at  257  Km. 
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Time  variation  of  the  6300  A  local-zenith  brightness  at 
100  and  200  km  from  the  exhaust  track's  center,  at 
331*  Km  at t i tude. 
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Figure  6a 
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Figure  6b.  Time  variation  of  the  6300  local-zenith  brightness  at 

100  and  200  km  from  the  exhaust  track's  center,  at  kOO  km 
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Figure  6c.  Time  variation  of  the  6300  A  local-zenith  brightness  at 

100  and  200  km  from  the  exhaust  track's  center,  at  1*50  km. 


Figure  6  shows  the  time  dependence  of  local-zenith  6300  A  radiance 
100  and  200  km  horizontally  from  the  exhaust  trajectory.  The  time  con¬ 
stants  for  decay  of  the  glow  are  longer  than  the  lifetime  of  the  0^D 
state,  which  indicates  (as  expected)  that  excitation  is  continuing  along 
these  sight  columns  during  the  data  period. 

Further  information  on  the  spatial  and  temporal  dependence  of  the 

i 

0  D  afterglow  in  the  first  M5  min  after  exhaust  of  and  H^O  in  the 
F  region  can  of  course  be  derived  from  Fig  b,  which  in  effect  presents 
all  the  radiance  information  present  in  the  sampling  of  low  light  level 
video  frames  recorded  from  the  aircraft. 


SECTION  2 


ATMOSPHERIC  SW I R  RADIANCE  STRUCTURE  MEASURED  FROM  AIRCRAFT  ALTITUDE 

INTRODUCTION 

During  April  and  September  1979,  and  August  1980,  the  Defense 
Nuclear  Agency  conducted  a  series  of  a i rcraf t -based  measurements  of  the 
space-time  structure  and  solar  illumination  dependence  of  yield  of  short- 
wavelength  infrared  rad i at i on  from  the  aurora  1  particle-excited  ionosphere. 
The  data  were  taken  from  AFGL's  Optical  Flying  Laboratory  (USAF  NKC-135A 
S/N  55-3120)  at  auroral  latitudes  near  70°  W  longitude  (1979,  flights 
from  Pea.e  AFB ,  NH)  and  1A7°  W  longitude  ( 1 58  0  ,  f  ron'  Eielson  AFB  ,  AK)  , 
with  a  cryocooled  radiometer  having  a  0.}6>°  (  6  millirad)  square  field 
of  viev.  pointed  near  the  magnetic  zenith.  The  principal  goal  was  to 
establish  the  spatial  and  temporal  correlations  of  2.8  -  }. 1. mi-band 
emission  with  the  deposition  of  ion i z i ng-par t i c 1 e  energy  in  the  dark  and 
sunlit  atmosphere,  at  a  1  Km  footprint. 

A  detailed  Test  Plan  for  the  aircraft  missions,  describing  the 
radiometer  with  its  coaligned  energy  input-referencing  air  fluorescence 
photometer  and  supporting  instruments,  and  outlining  a  procedure  for 
analyzing  the  sky  radiance  distribution  data,  is  in  Appendix  I  of  Ref  21. 
Measurements  of  the  atmosphere's  Rayleigh-scatter  background  when  particle 
deposition  altitudes  are  sunlit  are  presented  in  Section  5  of  Ref  1  , 
which  also  discusses  briefly  and  qualitatively  the  1979  flight  data.  A 
preliminary  review  of  the  data  set  was  given  in  Ref  22,  alone  with  back¬ 
ground  information  on  calibration  of  the  photometer. 

Summarized  results  from  earlier  flights  of  a  1 6 0 -m i 1 1 i rad-field 
SWIR  radiometer  and  photometer  (footprint  17  km  at  emission  altitudes) 
are  given  in  Ref  21.  In  this  Section  we  evaluate  the  higher  signal/ 
noise  data  segments  from  the  narrow-field  instruments,  focusing  on  the 


dependence  of  the  output-input  correlation  and  infrared  chemi 1 uminous 
yields  on  auroral,  sol ar- i 1 1 umi nat ion  ,  and  measurement  parameters. 

BACKGROUND 

The  source  simulating  excitation  by  radiations  from  nuclear  bursts 
is  the  energetic  (^5  keV)  electrons  from  the  magnetosphere  that  precipi¬ 
tate  into  the  h i gh- 1  a t i tude  upper  atmosphere  to  produce  the  optical  and 
infrared  aurora.  The  input  fluxes  and  peak  altitudes  of  energy  deposi¬ 
tion  are  measured  from  the  air  fluorescence  these  particles  excite  at 
visible  wavelengths.  Radiance  in  the  N*  First  negative  (0,0)  39 1 ^  A 
band,  which  is  known  to  be  closely  proportional  to  the  total  column 
energy  deposition,  is  monitored  by  a  photometer  whose  field  of  view 
coincides  with  that  of  the  2.832  -  3-125um  (FWHM) -sens i t i ve  radiometer. 
Spatial  distributions  of  input  in  the  vicinity  of  the  radiometer-photo¬ 
meter  fields  (and  auroral  type)  are  determined  each  1/30  sec  by  a  low- 
light-level  video  camera,  backed  up  by  a  photographic  all -sky  camera. 

A  second  radiometer  monitors  intensity  in  the  hydroxyl  vibrational 
(Meinel)  bands,  whose  fundamental  presents  a  very  s 1  owl y-va ry i ng  back¬ 
ground  in  the  SWIR  radiometer  signal,  which  is  believed  (Ref  21 )  to  be 
due  to  nitric  oxide  overtone  emission.  A  filter  photometer  measures 
ratios  of  emission  intensity  from  which  altitude  profiles  of  energy 
deposition  can  be  determined  (Appendix  II  of  Ref  23  ) ;  and  in  some  of 
the  missions  (in  1980)  altitude  profiles  of  ionospheric  electron  density, 
from  which  instantaneous  energy-input  profiles  can  also  be  estimated, 
were  measured  along  the  aircraft's  track  by  the  incoherent-scatter  radar 
at  Chatan ika  ,  AK. 

The  atmospheric  species  ionized,  dissociated,  and  raised  to  ex¬ 
cited  states  by  direct  primary  electron  impact  and  various  secondary 
processes  react  with  ambient  molecules  and  atoms,  producing  further 
excited  species  that  radiate  at  infrared  wavelengths.  Chemi luminous 
yields  ("efficiencies")  in  emission  features,  and  the  times  for  grow-in 
and  decay  of  the  infrared  radiations,  serve  (with  emission  spectra)  in 
identifying  the  radiating  species  and  the  exc i tat ion/de-exc i tat  ion 
processes  (a  further  goal  of  the  aircraft  program).  This  information! 

*»5 


is  needed  both  to  permit  scaling  of  simulation  data  to  the  altitude 
profiles  of  energy  deposition  that  result  from  atmospheric  nuclear 
explosions  and  to  determine  the  overlap  of  the  emission  spectrums  on 
wavelength  regions  to  which  surveillance  systems  are  designed  to  respond. 

Resolution  of  the  radiance  measurements  discussed  here  (£  1  km 
at  energy  deposition  altitudes,  and  1  sec  temporal)  is  intended  to  allow 
SWIR  output  and  "clutter"  to  be  correlated  with  energy  input  at  the 
scale  of  interest  to  spaceborne  sensors.  Quan t i f i ca t ion  of  the  visible- 
SWIR,  input-output  correlation  would  also  permit  use  of  the  existing 
extensive  photographic  data  base  on  visible  skyglows  resulting  from  past 
atmospheric  nuclear  tests  to  characterize  sky  background  variations  in 
the  2.8  -  3.1  urn  band. 

APPROACH 

The  atmosphere  above  the  subsonic  jet  aircraft's  operating  altitude 
of  '35  kft  is  essentially  transparent  to  downward-directed  radiation 
between  2.8  and  "3.1pm;  this  band  encompasses  about  half  of  the  emission 
in  the  NO  molecule's  overtone  (Av  =  2) sequence  as  excited  in  chemi- 
luminous  reactions  following  hard  particle  irradiation  of  air  (see,  for 
example,  Fig  2  of  Ref  1  ;  the  fractional  response  of  the  radiometer  to 
the  sequence  is  derived  in  Ref  21  ).  Initially,  the  atmospheric  emission 
at  these  wavelengths  was  measured  from  the  NKC-135A  aircraft  using  a 
radiometer  and  coaligned  photometer  having  1/6-radian  fields  of  view 
pointing  toward  the  zenith,  with  an  effective  signal  integration  time 
of  30  sec  (Ref  1  ) . 

Results  of  these  earlier  data  flights  (Ref  1  ),  which  refer  to  air 
between  MOO  and  150  km  altitude  in  which  M0^  ^  ion  pairs/(cm^  sec) 
are  being  produced,  can  be  summarized  as  follows. 

1.  SWIR  output  follows  particle-energy  input 
in  time  and  space  to  within  the  instru¬ 
ment's  resolution; 
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2.  The  output  of  SWIR  photons  increases  linearly 
with  input  at  auroral  dosing  rates; 

3.  Lumped-parameter  efficiency  for  emission  of 
what  is  interpreted  as  the  vibrational 
overtone  of  NO  shows  stat i st ical 1 y-s i gn i f i cant 
differences  among  both  long  (1/2  hr)  and  short 
(1  min)  segments  of  data,  with  yields  ranging 
from  0.42  to  0.72. 

These  results  from  the  wide-field  instrument  flights  are  by  and  large 
consistent  with  those  from  rocketborne  radiometers  pointed  toward  aurora 
(Ref's  1,  23),  and  (with  the  exception  of  the  as  yet-unexplained  variab¬ 
ility)  with  rate-coefficient  data  from  laboratory  investigations  of 
vibrational  excitation  of  the  nitric  oxide  molecule. 

These  earlier  flights  also  showed  that  measurements  of  both  the 
SWIR  emission  and  its  referencing  visible  air  fluorescence  have  adequate 
signal/noise  when  the  altitudes  at  which  precipitating  particles  deposit 
their  energy  are  sunlit  (as  reviewed  in  Section  5  of  Ref  1  ).  Operating 
near  the  equinoxes  and  navigating  to  compensate  for  wind  drift,  the 
aircraft  is  able  to  hold  the  solar  depression  within  the  requisite 
'^7 -10°  range  for  2-3  hours  while  staying  inside  the  (average,  Q  =  3) 
auroral  oval  and  permitted  airspace.  We  note  that  i dent i f i cat i on  of  day- 
night  differences  in  SWIR  dynamics  and  yield  would  indicate  a  deficiency 
in  current  understanding  of  the  processes  that  lead  to  vibrat ional ly 
excited  NO  in  part i c 1 e- i rrad i ated  air,  as  current  DoD  code  models  of 
the  sky  background  do  not  include  any  effects  of  excitation  by  solar 
photons.  (Additionally,  no  difference  from  night  to  twilight  (or  day) 
in  the  concentrations  of  ambient  species  that  affect  the  NC^ -produc i ng 
reaction  chain  can  be  identified.)  To  assess  the  behavior  of  SWIR 
emission  during  "sunlit"  aurora,  several  of  the  flights  of  the  narrow- 
field  instruments  were  made  at  the  twilight  transition,  with  auroral 
energy-deposition  altitudes  partially  or  totally  sunlit. 


The  narrow-field  instruments  are  installed  tilted  15°  forward  of 
the  static  airframe's  vertical,  so  that  at  auroral  latitudes  they  point 
in  the  direction  of  the  geomagnetic  field  when  the  fuselage  (which  is 
normally  pitched  2°  nose-up  in  flight)  is  moving  southward  along  a  mag¬ 
netic  meridian.  In  this  pointing  mode  the  precipitating  auroral  electrons 
within  the  fields  of  view  have  unique  initial  energy  spectrums  and  thus 
produce  a  defined  altitude  profile  of  atmospheric  excitation.  The  fields 
are  precisely  coaligned  (<  0.03°  estimated  error)  by  means  of  a  specially- 
designed  optical  bench,  construction  of  which  is  documented  in  Ref  24. 

An  independent  multi-channel  photometer  with  a  somewhat  wider  field  (2° 
full  angle)  is  installed  A3  feet  to  the  rear  of,  and  approximately  aligned 
with,  the  prime  instruments;  its  major  function  is  determination  of 
energy  deposition  altitudes  from  the  column  intensity  ratios  of  the  01 
6300  A  and  N+  4278  A  auroral  features  (as  described  in  Appendix  II  of 
Ref  23)<  (Calculations  of  the  collisional  quenching  of  N0  +  by  0  (Ref  25) 
in  fact  predict  a  minimum  of  SWIR  yield  near  100  km,  which  may  be  present 
in  some  of  the  rocket  radiometry  data  (Ref  1  )). 

ALTITUDE  PROFILE  DETERMINATIONS 

The  two  color  method  (Ref 23  and  earlier  references  therein)  quanti¬ 
fies  the  self-evident  principle  that  the  altitude  profiles  of  emission 
of  individual  molecular  bands  and  atomic  lines  are  functions  of  the 
atmosphere's  composition  profile,  and  so  reflect  the  penetration  profile 
and  thus  the  energy  distribution  of  the  primary  precipitating  auroral 
particles.  Major  assumptions  in  the  theory  are  1)  uniqueness  of  this 
energy  spectrum  within  the  photometer's  field  of  view,  2)  validity  of  its 
assumed  shape  (Maxwellian),  and  3)  steady  state  excitation/depopulation 
of  the  upper  state  of  the  0^D  -  ^P  6300  A  transition.  The  first  condi¬ 
tion  is  satisfied  during  segments  of  flights  in  which  the  aircraft  fusel¬ 
age's  heading  is  near  180°  (magnetic);  but  in  others,  particularly  the 
twi 1 i ght-trans i t ion  missions  where  the  trajectory  is  chosen  to  maintain 
constant  solar  depression,  the  angle  between  the  geomagnetic  zenith  and  the 
instrument's  field  can  be  as  high  as  26°  (13°  +  13°;  MOO  km  horizontal 
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separation  between  the  and  01  red  line  emission  altitudes). 

Thus  the  confidence  placed  on  altitudes  derived  by  this  method  depends  on 
both  aircraft  heading  and  the  scale  length  of  uniformity  of  the  aurora- 
producing  "beam."  The  time  to  reach  steady-state  conditions  in  01  6300  A 
is  at  least  one  lifetime  of  the  upper  state  against  radiative  and 
collisional  depopulation,  which  is  about  30  sec;  thus  the  method  provides 
more  reliable  altitudes  of  auroral  forms  that  are  spatially  and  temporally 
constant  at  this  scale. 

Additionally,  deviations  from  the  assumed  shape  of  the  particle 
energy  spectrum  may  be  a  potential  source  of  error  in  the  deposition 
altitude  prof i 1 e .  1  n  consequence  the  1980  flights  were  coordinated  with 
Chatanika  radar  scans  which  determined  altitude  profiles  of  ionospheric 
electron  densities,  from  which  input  rates  can  be  derived  with 
the  assumption  of  constant  deposition  rate  over  the  previous  M0  sec 
(refer  to  Ref's  1  and  23  for  a  review  of  the  measurement  concept).  A 
total  of  some  15  hr  of  coordinated  measurements  in  5  data  missions  was 
made,  with  the  aircraft  flown  back  and  forth  along  (or  near)  the  radar's 
magnetic  meridian.  An  attempt  was  made  to  synchronize  the  tracking 
radar's  elevation  scan  with  the  aircraft's  trajectory  so  that  the  radio¬ 
meter  and  photometer  fields  would  overlap  the  10  millirad  radar  field  at 
100  km  altitude.  Flight  profiles  that  would  obtain  an  adequate  sampling 
of  both  sunlit  and  nighttime  aurora  were  designed  for  each  mission,  with 
the  aircraft  entering  the  data  taking  track  when  the  solar  depression 
reached  7°  at  evening  twilight. 

The  instruments'  100-km  intercept  was  flown  1°  02'  A3"  N  and  S 
along  the  meridian  through  Chatan i ka  (between  64°  03'  31"  and  66°  08'  53"), 
which  is  a  total  distance  of  266  km  when  the  29°  magnetic  declination  is 
considered.  At  the  extremes  of  this  track  the  radar  beam's  zenith 
angle  is  49. 4°,  and  its  angular  speed  is  0.056  degrees/sec  (it  is  0.13 
deg/sec  when  the  a i rcraf t  i  s  in  the  zen i th)  .  The  signal  integration  time 
is  a  factor  }-£>  greater  than  in  a  previous  a i rc raf t- rada r  coordination 
(Ref  1  ).  The  number  of  two-way  meridian  passes  varies  between  3  and  6 
i n  the  1 980  f 1 i ghts . 
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Since  operational  considerations  required  that  the  radar  facility 
be  given  a  preprogrammed  tracking  routine, no  provision  could  oe  made 
in  advance  for  the  effect  of  winds  on  the  aircraft's  position.  There¬ 
fore  some  offsets  developed  between  the  two  fields.  ( Further,  c rosswi nds 
frequently  result  in  "crabbing"  of  the  airframe,  which  offsets  the  in¬ 
strument  fields  from  the  predetermined  flight  track.)  To  mitigate  this 
problem  a  series  of  checkpoints  at  wh:ch  the  aircraft  would  arrive  was 
set  up.  The  aircraft  was  not  permitr''d  to  deviate  from  its  course  along 
the  predetermined  track  to  "chase"  auroral  forms,  which  of  course  lessened 
the  probability  of  its  encountering  intense  particle  bombardment. 

Electron  densities  along  the  scanning  beam  have  been  received 
from  the  radar  group;  these  have  not  yet  been  converted  to  the  contour 
plots  (of  the  type  shown  in  Fig's  20  and  21  of  Ref  1  )  needed  to  deter¬ 
mine  the  energy  input  rates  within  the  aircraft  i ns t rumen t s ' f i e 1 d  of 
view.  We  note  that  the  inherent  ambiguity  (time-space  mix)  of  the 
I  energy  input  rate  profiles  derived  from  elevation  scans  from  a  single 

instrument  location  limits  the  application  of  t'e  radar-backscatter 
data  in  interpreting  the  SWIR  radiation  data,  as  the  spatial  distribu¬ 
tion  of  input  flux  can  change  before  the  radar  beam  completes  its 
traverse  across  the  column  defined  by  the  instrument  field  (refer  to 
the  comments  in  the  report  on  the  earlier  coordination,  Ref  1  ).  Pend¬ 
ing  reduction  of  these  aircraft  data  and  an  assessment  of  effective 
offset  of  fields,  the  peak  energy  deposition  altitudes  presented  here 
have  been  ierived  using  the  spectroscopic-ratios  method. 

DATA  REDUCTION 

A  total  of  15  high  signal/noise  data  segments  from  10  fl ights .which 
included  70  min  of  measurements,  was  selected  from  the  data  base  for 
the  initial  analysis.  Auroral  precipitation  intensity  ("input")  was 
typically  >  IBCII+,  and  the  segments  were  selected  to  include  totally 
sunlit,  partially  sunlit,  and  totally-dark  periods.  A  typical  set  of 
calibrated  radiance  traces,  showing  the  391^*  A  band  fluorescence, 

2.8  -  3.1pm,  and  Av  =  2  hydroxyl  is  in  Fig  7  (01  5577  A  intensities 
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measured  by  the  narrow-angle  photometer  are  not  shown).  The  wider-field 
OH  overtone  (1.66  -  1.7^um  FWHM)  a i rg low  monitor's  optic  axis  pointed  to 
the  zenith  in  1979  and  was  nominally  coaligned  with  the  prime  instru¬ 
ments  in  some  of  the  1 980  flights.  The  ratio  of  in-band  OH  fundamental 
to  measured  overtone  emission  intensity  has  been  calculated  to  be  about 
2.5  :  1;  in  practice,  variations  in  OH  background  over  the  short  data 
segments  are  so  small  that  they  do  not  affect  the  analysis. 

The  fluorescence  photometer  installed  in  1980  was  a  new  design 
whose  initial  calibration  is  only  indirectly  applicable  to  the  measure¬ 
ments  made  in  1979.  A  final  calibration  for  1979  data  was  made  by 
cross-correlating  the  raw  output  voltage  readings  from  the  instrument 
with  the  absolute  intensities  that  were  measured  by  the  aircraft's 
12-channel  filter  photometer  (Ref  1  ) ,  as  described  in  Section  k  of 
Ref  22  .  No  instrument  performance  irregularities  were  encountered  in 
the  1 980  mi ss ions . 

Measurement  and  geophysical  conditions  for  the  15  traces  selected 
are  shown  in  Table  k  ,  and  montages  of  all -sky  photographs  (165°  circular 
field,  sec  on  EK2f* 75  at  f/2)  of  the  auroral  input  distribution  are 
in  Fig  8  .  The  negatives  for  each  segment  have  all  been  printed  with 
the  same  exposure  to  indicate  relative  scene  brightnesses.  Segment  10 
has  been  omitted  from  Fig  8  because  the  signal/noise  ratio  in  its 
radiance  traces  was  found  to  be  below  that  acceptable  for  analysis. 

The  entries  in  Tableware  as  follows.  The  first  column  lists 
the  identifying  numbers  assigned  to  each  data  segment.  Next  is  flight 
date  and  time  of  meridian  passage  at  the  midpoint  of  the  data  segment 
(which  is  needed  to  determine  solar  elevation),  whose  duration  is  in 
the  third  column.  The  next  three  columns  give  aircraft  position, 
fuselage  heading,  and  direction  of  the  flight  trajectory;  the  latter 
two  azimuths  typically  differ  by  a  few  degrees  because  of  crosswinds. 

(When  the  fuselage  is  oriented  at  180°  magnetic,  the  coaligned  input- 
and  output-measuring  instruments  point  directly  up  the  field  lines; 
yaw  of  the  airframe  by  azimuth  angle  4  results  in  a  misalignment  of 
only  4  sin  13°  =  0.2254'.) 
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Table  4.  Auroral  and  observation  parameters  during  SWIR  enhancements 
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Figure  8c.  All-sky  montages  of  the  auroral  energy  input  distribution 
before  and  during  data  segment  3. 
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Figure  8g.  All-sky  montages  of  the  auroral 
before  and  during  data  segment  7 
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Figure  8h.  All-sky  montages  of  the  auroral  energy  input  distribution 
before  and  during  data  segment  8. 
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Figure  8j .  All-sky  montages  of  the  auroral  energy  input  distribution 
before  and  during  data  segment  11. 
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All-sky  montages  of  the  auroral  energy  input  distribution 
before  and  during  data  segment  12. 
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Figure  8n.  All-sky  montages  of  the  auroral  energy  input  distribution 
before  and  during  data  segment  15  (1  of  3  pages). 
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-igure  8n  (concluded).  All-sky  montages  of  the  auroral  energy  input 

distribution  before  and  during  data  segment  15  (3  of  3  pages). 
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The  elevation  angle  of  the  center  of  the  solar  disk  at  the  aircraft 
and  at  the  latitude  and  longitude  of  the  120-km  altitude  instrument 
intercept  the  directly  applicable  figure  --  follows  the  solar  declina¬ 
tion  (which  has  been  tabulated  because  it  is  an  input  to  the  program 
for  calculating  solar  elevation).  Symbols  D  and  I  indicate  decreasing 
(evening  twilight)  or  increasing  (morning  twilight)  elevation,  to 
indicate  whether  the  region  of  atmosphere  on  which  the  optical  measure¬ 
ments  are  made  had  been  sunl i ght-"predosed ."  The  column  following  is 
the  hard  earth's  shadow  height  at  the  13°  zenith  angle  intercept  (as 
defined  in  Fig  27  of  Ref  1  ).  The  tangent  altitude  at  which  the  ( 1 966 
U.S.  Standa rd) atmosphere  transmits  half  the  incident  solar  ultraviolet 
radiation  is  34  km  at  3400  A,  45  km  at  3050  A,  and  70  km  at  2380  A; 
the  minimum  auroral  altitude  reached  by  half  the  photons  at  these  UV 
wavelengths  can  be  found  with  adequate  accuracy  by  adding  these  "screen¬ 
ing"  heights  to  the  illuminated  heights  listed  in  Table  4  . 

Column  intensity  ratios  of  visible-aurora  features  at  the  times 
in  the  data  segment  when  energy  input  rate  is  a  maximum,  are  listed  in 
the  next  column.  These  are  followed  by  the  altitudes  of  peak  energy  de¬ 
position  at  these  more  or  less  characteristic  times,  calculated  from  the 
6300  A/4278  A  ratio  and  4278  A  intensity.  These  altitudes  refer  to  the 
aircraft's  zenith  in  Flights  908  anJ  909  only,  after  which  the  field 
of  view  of  the  wave  1  eng th-cyc 1 i ng  photometer  was  directed  13°  forward. 

The  final  two  columns  list  the  three-hour  local  (K)  and  global 
(Kp)  magnetic  indexes  at  the  times  of  the  flight  segments.  Local  K's 
refer  to  the  College,  AK  sector  and  thus  may  not  be  directly  applicable 
to  the  1979  flights,  which  were  80°  in  longitude  east  of  College. 

The  fractional  solar  illumination  of  each  energy-deposition 
column,  which  is  listed  in  Table  5,  was  determined  as  follows.  We 
first  assigned  an  "input"  altitude  profile  based  on  the  peak-deposition 
altitudes,  following  the  procedure  in  Ref  26.  This  altitude  distribution 
was  taken  as  remaining  fixed  for  the  typically  2-min  duration  of  the 
SW I R-measurement  interval.  (For  Segments  15a  and  b,  an  average  peak 
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Table  5 


Fractional  illumination  of  the  altitude  profile  of  energy 
deposition  in  the  15  data  segments.  Values  refer  to  the 
"hard  earth"  definition  of  solar  illuminated  height,  as 
given  in  Fig  27of  Ref  1. 


%  SUNLIT  . -  HARD  EARTH  DEFINITION 
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altitude  somewhat  higher  than  that  listed  in  Table  5  was  used.)  We  then 
estimated  the  fraction  of  the  energy  deposition  lying  above  the  shadow 
height  from  ratio  of  areas  under  the  profile  curve.  Six  segments  are 
identified  as  being  essentially  totally  sunlit  (3  each  from  the  1979 
and  1980  series),  six  further  are  partially  illuminated  (again  3  from 
each  year),  and  the  remaining  three  are  nighttime  aurora.  Reference  is 
made  to  the  fact  that  the  entries  in  Table  5  do  not  consider  outscatter ing 
of  solar  photons  by  the  lower  atmosphere;  as  noted,  scattering  results 
in  lower  fractional  illumination  by  ultraviolet  photons. 

DATA  ANALYSIS 

The  aircraft  traverses  the  projection  of  the  instrument  fields  at 
the  altitudes  where  incoming  particles  deposit  energy  in  about  4  sec, 
and  the  integration  time  needed  to  achieve  adequate  SWIR  signal /noise 
is  somewhat  less.  If  the  time  required  for  aerochemical  reactions  to 
excite  SWIR  radiation  is  less  than  or  comparable  to  4  sec,  or  under 
conditions  that  production  and  loss  rates  have  reached  equilibrium  in  the 
atmospheric  volume  probed,  the  cross-plot  method  that  we  applied  to  the 
wide-field  data  (Ref  1  )  is  a  valid  means  for  determining  average  SWIR- 
chemi luminous  yields.  We  therefore  made  cross-plots  of  the  simultaneously 
measured  "input"  and  "output"  radiances,  reading  the  traces  at  5_sec 
intervals  in  cases  where  the  duration  of  enhancement  was  several  min  and 
i  to  1  sec  when  it  was  M  min  (Segment  12,  for  example).  Limited  manual 
smoothing  was  applied,  to  minimize  bias  against  short-period  fluctuations 
while  ensuring  that  noise  spikes  are  suppressed.  As  the  contribution  to 
the  SWIR  signal  from  OH  airglow  remains  essentially  constant  over  indi¬ 
vidual  data  segments,  it  was  not  subtracted;  this  of  course  does  not 
change  the  slope  of  the  cross-plots. 

Example  cross-plots  from  three  data  segments  representing  one  day 
(Segment  *12)  ,  one  night  (Segment  1)  and  one  twilight  (Segment  3)  I BC  11  + 
aurora  are  in  Fig  9  .  The  radiance  traces  had  been  filtered  in  the 
previous  computer  merging/calibration  of  the  original  samplings,  each 
plotted  point  being  an  average  of  11  points  spaced  0. 01  apart.  For 
Segment  12  our  sampling  interval  was  0.4  sec,  and  for  1  and  3,  1.0  sec. 
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The  visible-auroral  distributions  in  all-sky  photographs  can  best  be  de¬ 
scribed  as  a  near-uniform  arc  in  Segments  3  and  12,  with  Segment  1  a 
non-stationary  "breakup." 

As  reported  previously  (Ref  1)  some  of  the  SWIR  enhancements 
show  several-sec  longer  decay  times  or  offsets  of  emission  peaks  from  the 
39 1 ^  A  fluorescence  signal  (Segment  12  is  an  extreme  example).  The  "loop¬ 
ing"  of  the  cross-plots  in  Fig  9  is  one  result  of  the  decorrelat ion  of  the 
two  column  radiances;  for  example  in  segments  where  the  energy  input  is 
decreasing  the  "delayed"  i ns tantaneous  SWIR  radiances  lie  above  the  mean 
regression  line,  as  is  particularly  evident  in  Fig  9b  and  c.  This 
effect  appears  in  both  sunlit  and  nighttime  flights,  as  discussed  further 
below.  We  have  "averaged"  the  data,  also  as  discussed  later,  in  deriving 
mean  chemi 1 umi nous  yields. 

The  straight  lines  in  Fig  9  have  been  least-squares  fitted  to  the  data 
points,  using  the  procedure  and  criteria  applied  to  the  wide-field  instru¬ 
ment  data  in  Ref  21.  The  R-values  represent  statistical  correlation  co¬ 
efficients,  or  "goodness  of  fit"  with  a  straight  line,  and  error  given  is 
the  30%  confidence  limit  of  the  derived  1  ine  slope.  We  have  again  assumed 
that  the  uncertainty  in  the  fluorescence  intensities  in  the  391^  A-band  is 
small  compared  to  that  of  the  SWIR  band  signal.  The  assumption  of  linearity 
of  course  does  not  hold  when  the  fall'  time  of  SWIR  radiance  exceeds  that  of 
391A  A,  as  in  Fig  9c;  note  that  the  seven  circled  points,  which  refer  to 
the  first  2i  sec  of  the  8  sec  measurement  interval  when  the  instrument 
field  first  entered  a  stable  arc  and  the  radiances  are  increasing,  do  show 
a  high  linear  correlation  when  considered  separately. 

Under  the  conditions  stated  above,  the  slope  of  the  "best  fit" 
straight  line  can  be  interpreted  as  being  proportional  to  the  chemilumin- 
ous  energy  efficiency  (or  yield).  In  cases  such  as  Fig  9c  the  fit  to  all 
the  data  points  makes  a  reasonably  accurate  estimate  of  the  average  SWIR 
yield  over  the  measurement  interval  since  about  an  equal  number  of  points 
lie  above  and  below  the  regression  line  (the  correlation  coefficient  pro¬ 
vides  a  measure  of  this  accuracy).  This  assumption  is  further  supported 
by  results  of  calculations  of  the  ratio  of  "area"  under  each  SWIR  enhance¬ 
ment  (in  kR-sec)to  that  under  its  391 4  A  radiance  increase,  which  show 
the  time-integrated  photon  emission  rate  ratios  to  be  in  agreement  with 


the  best-fit  slopes  (within  the  limits  of  the  precision  of  the  data).  This 
second  method  for  determining  yield  requires  subtraction  of  the  OH-funda- 
mental  airglow  from  the  SWIR  radiances.  This  airglow  intensity  was  assessed 
by  interpolation  of  the  background  on  either  side  of  (before  and  after)  the 
aurora-associated  enhancement.  The  uncertainties  in  slope  listed  on  Fig  9 
refer  only  to  90%  statistical  confidence  limits,  and  do  not  reflect  other 
potential  systematic  and  erratic  sources  of  error  due  to  calibrations, 
instrument  instability,  window  frost,  and  similar  factors. 

Following  the  procedure  that  we  applied  in  Ref  21,  and  taking  0.42 
as  the  fraction  of  NO  Av  =  2  sequence  within  the  radiometer's  bandpass 

(usinq  the  COCHISE  population  model),  we  determined  chemi luminous  yields 
from  the  slopes  of  each  of  the  15  segments'  cross-plots.  (in  three  cases 
the  "area"  method  was  used,  as  identified.)  The  results,  listed  in  order 
of  increasing  solar  depression  angle,  are  shown  in  Table  6  . 

The  column  labeled  PhM  delay  in  Table  6  gives  a  semi -quant i tat i ve 
measure  of  the  delay  between  SWIR  output  and  energy  input,  which  we  deter¬ 
mined  by  visually  comparing  offsets  of  emission  peaks  and  lags  in  decay 
time  between  the  two  data  traces.  The  estimation  method  was  checked  by 
applying  it  to  the  5577  A  (01  ' S  -  *0)  and  3914  A  photometer  signals;  the 
lag  of  the  metastable  oxygen  emission  was  less  than  1  sec  in  all  cases, 
which  is  larger  than  its  expected  value  of  sec  (Ref  13).  A  mean  char¬ 
acteristic  filtering  time,  which  refers  to  both  the  rise  and  fall  of  the 
SWIR  signal  and  which  is  derived  (by  M.  Bruce  of  AFGL)  from  Fourier  trans¬ 
forms  of  the  individual  data  streams,  is  given  under  the  heading  MB  LAG. 
Note  that  this  time  is  not  the  lag  derived  from  cross-correlation  of  the 
two  data  sets,  discussed  in  Appendix  I  of  Ref  21;  rather,  it  is  the  width 
of  the  low-pass  temporal  filter  that  "fits"  the  SWIR  output  to  the  input's 
frequency-space  components.  The  final  evaluation  of  the  aircraft  data  will 
apply  cross-correlation  and  delay-weighting  methods  to  quantify  the  induc¬ 
tion  time  of  chemiluminescence. 

Data  quality  (based  on  signal/noise  of  the  SWIR  radiance)  is  given 
in  the  next  column;  the  times  in  parentheses  are  the  lengths  of  the  en¬ 
hancement  segments.  Fuselage  pointing  and  altitude  of  peak  energy  deposi¬ 
tion  are  then  repeated  from  Table  4  to  permit  a  rapid  comparison  between 
yield  and  lag  and  these  measurements  parameters.  "Averages"  shown  for 
some  segments  refer  to  peak  emission  altitude  averaged  over  the  measure- 
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Table  6.  Chemi I uminous  yields  and  apparent  SWIR  delays,  with  deposi¬ 
tion  altitude  and  aircraft  heading  from  Table  4.  The  data 


Averaged  over  entire  data  segment. 
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1  i  sted 
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merit  period;  note  that  these  are  higher  than  the  altitudes  at  maximum  en¬ 
hancement,  as  the  sprectral  index  (or  hardness)  of  the  precipitating  par¬ 
ticles  tends  to  decrease  as  their  total  flux  decreases. 

The  efficiencies  listed  in  Table  6  show  a  statistically  signifi¬ 
cant  variability  of  more  than  a  factor  2  over  the  18  month  measurement 
period,  among  individual  flights  of  the  three  series,  and  even  within 
single  flights  (for  example , Segments  14  and  15a-b,  Flight  027),  A  com¬ 
parably  large  variation  was  present  in  the  yields  measured  by  the  wide- 
field  instruments  in  1975“78  (Ref  1  ),  whose  average  was  found  to  be 
about  i%.  With  the  narrow-field  instruments  the  mean  energy  yield  in 
SW| R-band rad i at  ion  in  part i cl e- i rrad iated ,  auroral  altitude  air  is  closer 
to  1%,and  in  some  cases  this  yield  reaches  1.A$, 

The  "lag"  between  SWIR  radiation  output  and  energy  input,  as 
estimated  from  offsets  between  emission  peaks  and  delays  in  the  SWIR 
signal's  return  to  baseline  after  the  excitation  has  essentially  dis¬ 
appeared,  also  varies,  showing  values  between  zero  and  6  sec.  (The 
Fourier-components  rat io  calcul at  ion  shows  filtering  time  constants  that 
cluster  closely  around  1  sec.)  It  is  instructive  to  note  that  no 
negative  lags  are  identified  by  either  method. 

CORRELATION  WITH  AURORAL  AND  MEASUREMENT  PARAMETERS 

To  identify  probable  sources  of  this  variability  in  both  mean 
yield  and  lag  time  of  SWIR  emission,  we  c i ^ss-plotted  results  against 
auroral  and  measurement  parameters  (Fig's  10  and  11).  The  circled 
numbers  in  the  plots  identify  the  data  segments. 

FiglOashows  no  apparent  correlation  of  yield  with  incoming  flux, 
which  is  not  unexpected  because  the  original  input-output  cross-plots 
are  closely  linear.  Results  from  the  earlier  wide-field  instruments 
support^this  finding..  Figure  1  Ob  indicates  that  SWIR  yield  is  also  in¬ 
dependent  of  fractional  auroral  illumination;  i.e.,  no  day-night  effect 
is  found.  No  correlation  is  found  with  global  magnetic  activity  index, 
or  altitude  of  the  aircraft  (not  shown),  on  which  could  depend  effects 
of  window  frosting,  atmospheric  attenuation  and  cirrus  (or  unexpected 
aspects  of  instrument  performance). 
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AVERAGE  NO  OVERTONE  YIELD,  % 

Figure  10e-f.  Average  NO  overtone  yield  plotted  against  altitude 
of  energy  deposition,  (e)  All  segments,  at  maxi¬ 
mum  depositon  rate;  (f)  south-pointing  segments, 
averaged  over  deposition  rates. 


fikft 


Some  evidence  for  a  trend  with  the  ratio  of  intensities  in  the 
hydroxyl  fundamental  and  overtone  nightglow  is  in  Fig  10c;  however  the 
negative  slope  (-3.23  ±  2.2  Av  =  1/Av  =  2  (in  photons)per  \%  change  in  ab¬ 
solute  SWIR  yield)  is  opposite  to  what  might  reasonably  be  expected.  Variable 
sensitivity  of  the  SWIR  radiometer  would  drive  boththis  hydroxy  1 -bands 
ratio  and  the  apparent  aurora-assoc iated  yield  in  the  same  direction; 
and  cirrus  clouds  above  the  aircraft  increase  the  apparent  ratio  (excess 
thermal  radiation  being  mistaken  for  airglow)  without  greatly  impacting 
the  yield  (i  ncreased  outscatter i ng  of  shorter  wavelength  photons  compen¬ 
sating  the  increased  attenuation  of  SWIR  radiation  from  the  aurora).  In 
any  case,  Fig  10c  has  too  low  a  correlation  coefficient  and  statistical 
precision  for  a  compelling  argument  supporting  the  ant i -cor rel at i on 
(note  that  the  trend  line  is  significantly  "weighted"  by  two  data  points, 
from  Segments  12  and  14). 

Fig  lOe  may  be  interpreted  as  suggesting  a  correlation  between 
average  SWIR  yie'd  and  the  altitude  of  peak  energy  deposition  by  the 
incoming  electrons.  The  altitudes  in  Fig  lOe  refer  only  to  that  section 
of  each  data  segment  at  which  the  39 1 2f  A  column  brightness  (i.e.,  in¬ 
coming  particle  flux)  was  a  maximum.  For  the  five  cases  in  which  the 
aircraft  was  heading  magnetic  south  (Fig  lOf),  peak  energy  deposition 
altitudes  were  averaged  over  the  full  measurement  interval.  The  slope 
of  the  regression  line  fit  to  all  data  points  is  3.6,  which  would  be 
an  increase  in  NO  overtone  efficiency  of  about  0.28%/km  altitude;  the 
low  correlation  coefficient,  0.08,  indicates  a  very  poor  linear  fit, 
and  furthermore  the  90%  confidence  limits  on  the  slope  are  very  broad. 

If  Segment  7  is  omitted,  the  slope  increases  to  16.8  (0.06%/km  altitude) 

R  =0.49.  For  both  sample  sets,  the  large  statistical  uncertainty 
indicates  that  the  slope  of  the  best  fit  line  could  result  in  SWIR 
efficiency  increases  as  low  as  ±0. 03%/km  of  peak  energy  deposition 
al t i tude . 

To  investigate  further  this  apparent  increase  in  SWIR  yield  with 
deposition  altitude,  "instantaneous"  efficiencies  and  altitudes  were 
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Table  7.  Correlation  of  SWIR  enhancement  delay  with 
auroral  and  measurement  parameters 


SOLAR  DEPRESSION 

<  10° 

10°  to  12 
12°  to  14° 

>  14° 

IBC  CATEGORY  OF  AURORA 

I  +  I  I 
III  +  IV 

I  +  I  I  +  I  I  I 
IV 

A/C  POINTING 

N 

S 

SWIR  YIELD 

>  n 

<  n 

PEAK  DEPOSITION  ALTITUDE 

<  120  km 
120  km 

>  120  km 

AURORAL  FORM 

Isolated  Arc  (IA) 
Other  (0) 

IA  North  Heading 
IA  South  Heading 
0  North 
0  South 

OH  (A v  »  1)  RADIANCE 

<  300  kR 
>  300  kR 

K£ 

1  -  4 

5  -  8 


--  PhM  DELAY-- 
4  2  sec  >  2  sec 


1 

2 

1 

2 


3 

2 

2 

1 


3 

3 

4 
2 


6 

2 

7 

1 


1 

4 


3 

5 


5 

2 


2 

3 


5 
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3 
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3 
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measured  at  5  sec  intervals  for  the  3  k  mi n  of  Segment  13,  during  which 
a  temporal  1 y-un i form  (sunlit)  auroral  arc  was  traversed  by  instrument 
fields  pointing  in  the  direction  of  the  geomagnetic  field  (Fig  8).  This 
segment  represents  one  type  of  optimum  data:  an  isolated,  stable  (within 
the  6300  A  line  equilibrium  criterion)  arc  viewed  long  the  directions  of 
incoming  particle  precipitation.  The  results,  which  generally  support 
those  in  Fig  lOe, suggest  an  increase  in  mean  NO  overtone  yield  of  about  0.15& 
per  km  a  1 1 i tude, wi th  its  statistical  limits  0.33  and  0.061  per  km.  This 
finding  of  a  large  variation  with  altitude  does  not  agree  with  calcula¬ 
tions  shown  in  Fig  12  based  on  a  simple  model  of  N0^  exc i tat i on-deexci ta- 

2 

tion:  quenching  of  N0(v)  by  0  proportional  to  v,  quenching  of  N  D  by  0 

“13  3 

with  rate  coefficient  5  x  10  Jcm /sec,  [0]  profile  as  measured  from 

2 

DNA/AFGL  rocket  I R  807 . 57~  1,1 i  N  D/ion  pair  at  all  altitudes  (the  model 
is  based  generally  on  ideas  in  Ref  25 )  .  These  simplified  calculations 
indicate  that  the  vertical  column  integrated  NO  overtone  yields  should 
increase  very  slowly  with  peak  energy  deposition  altitude,  roughly 
0.0021  per  km  over  the  100-140  km  altitude  range. 

We  have  separated  the  estimated  SWIR  lags  into  two  groups,  "short" 

(<:  2  sec)  and  "long"  (>  2  sec),  for  comparison  with  some  of  those  auroral 
and  measurement  parameters  perceived  as  potentially  impacting  the  data. 

Table  7  gives  numbers  of  segments  falling  into  the  categories 

--  solar  depression  angle  (sunlit  aurora,  partially 
sunl i t ,  dark)  ; 

--  IBC  intensity  of  the  aurora; 

--  direction  in  which  the  radiometer  is  pointing 
fciarallel  to  or  at  26°  to  the  magnetic  field 
1 ines) ; 

--  SWIR  yield; 

--  altitude  of  maximum  energy  deposition; 

--  type  of  auroral  form,  as  determined  from  the 
a  1 1 -sky  photographs ; 


--  form  +  instrument  pointing; 

--  mission  dates; 

--  zenith  radiance  in  the  segment  of  hydroxyl 
fundamental  within  the  radiometer's  wave¬ 
length  sensitivity  band;  and 

--  global  magnetic  activity  index  Kp. 

It  is  apparent  that  there  is  no  statistically  significant  favoring  of 
long  or  short  delays  under  any  of  these  i nput-cond i t i ons  headings,  as 
within  the  expected  statistical  fluctuation  an  equal  number  of  cases  is 
included  in  each  category. 

Note  in  particular  the  lack  of  definite  correlation  of  lag  with 
solar  depression  angle  or  (more  importantly)  fractional  illumination  of 
aurora.  Fig  11,  which  cross-plots  the  delays  against  illumination  (and 
also  peak  excitation-rate  altitude),  illustrates  this  finding:  no  trend 
is  apparent.  Putting  the  data  samples  in  the  quantized  format  of 
Table  7 ,  we  have 


%  Sun  1  i  t 
100% 

<  100%,  >  0% 
0% 


--  PhM  Delay 

N<  2  sec  >  2  sec 

2  4 

3  3 

2  1 


Hard  earth  shadow  height  definition;  see  Table  5- 

We  may  conclude  that  the  data  sample  is  too  small  to  support  the  hypo¬ 
theses  of  a  variation  in  time  delay  of  SWIR  radiation  (as  qualitatively 
determined)  with  fractional  solar  illumination  or  mean  energy  deposition 
altitude.  Similarly,  the  mean  yield  of  SWIR  photons  in  a  particle- 
bombardment  event  does  not  correlate  with  solar  illumination,  or  con¬ 
versely,  a  hypothesis  of  a  day-night  difference  in  chemi 1 umi nous 
efficiency  is  not  supported  by  this  set  of  a i rcraft- instrument  data. 
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Note  that  the  set  of  parameters  in  Table  7  does  not  address  the 
issue  of  whether  the  change  in  energy  input  rate  within  the  photometer's 
field  is  dominated  by  the  aircraft's  motion  across  temporal  1 y-spat ia 1 1 y 
stable  particle  bombardment  or  by  variability  of  this  bombardment.  In 
stationary  auroral  forms  the  SWIR  lag  would  be  expected  to  be  lower,  as  it 
depends  only  on  transport  (diffusion  and  wind  drift)  of  the  intermediate 
reacting  species,  as  discussed  in  the  Test  Plan  (Ref  21);  This  could  in 
fact  make  tie  SWIR  radiance  appear  to  lead  the  energy  input.  The  video 
images  of  the  auroral  distribution  can  be  used  to  answer  the  question  of 
whether  the  lag  is  dominated  by  chemical  or  transport  effects. 

SUMMARY  AND  RECOMMENDATIONS 

The  set  of  measurements  from  subsonic  jet  aircraft  altitude  with  a 
sensitive  narrow-field  radiometer  and  particle  energy  i nput-referenc i ng 
photometer  is  a  prime  source  of  information  on  the  spat ial -temporal 

correlation  between  SWIR  emission  output  and  ion i z i ng-exc i t ing  input  < 

into  the  atmosphere  above  MOO  km.  The  auroral  2.8  -  3.1pm  radiances 

in  near-zenith  projections  measured  from  AFGL's  Flying  Optical  Laboratory 

over  a  period  of  hours,  complement  those  in  earth-limb  projections  to 

be  measured  over  a  period  of  minutes  from  DNA/AFGL's  forthcoming  ELIAS 

sounding  rocket.  Additionally,  the  only  experimental  verification  that 

the  equilibrium  intensity  and  grow-in  decay  characteristics  of  this 

radiation  are  the  same  when  the  deposition  region  is  sunlit  or  dark 

comes  from  the  aircraft  measurements. 

The  principal  findings  of  this  evaluation  of  the  "best"  segments 
of  the  1 97 9~ 1 980  aircraft  data  set  are  that  the  SWIR  radiance  follows 
column  energy  input  to  within  a  few  sec  (the  quantitative  temporal/ 
spatial  dependence  has  not  yet  been  extracted  from  the  data)  ,  and  that 
the  mean  energy  efficiency  of  excitation  of  SWIR  emission  varies  between 
about  0.6%  and  1.3%  with  no  dependence  on  solar  illumination  so  far 
detectable.  15  such  segments,  each  representing  typically  2  min  of 
measurements  on  IBC  II  +  aurora  of  various  forms,  were  evaluated 
systematically.  The  "equilibrium"  chemi 1 uminous  SWIR  yield  shows 
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stat  i  st  ical  1  y-s  i gn i  f  i  cant  var  i  ab i  1  i  tyi  wi  th  a  mean  of  (assuming  the 
auroral  radiation  enhancements  to  be  the  NO  overtone,  with  a  spectrum 
as  determined  from  AFGL  laboratory  and  HIRIS  data).  This  is  almost  a 
factor  2  higher  than  the  average  determined  f  om  aircraft  measurements 
with  a  coaligned  radiometer  and  photometer  having  a  much  wider  field. 

Both  sets  of  data  show  factor-two  variability  in  yield  even  among  differ¬ 
ent  periods  of  individual  flights  in  the  auroral  oval  (as  well  as  among 
flight  series),  which  considering  also  the  lack  of  correlation 
with  the  near-dc  zenith  sky  background  from  OH  v i brat ional -fundamental 
radiation,  makes  it  unlikely  that  this  variability  is  an  artifact  of 
the  instrument's  response  or  laboratory  calibrations. 

The  data  provide  some  evidence  suggesting  that  SWIR  energy  effic¬ 
iencies  integrated  over  the  near-vertical  measurement  column  depend  on 
the  altitude  profile  of  energy  deposition  by  the  precipitating  auroral 
electrons  (Fig's  lOe  and  12).  While  such  an  increase  with  peak  disposition 
altitude  might  be  expected  from  the  increased  quenching  of  NO^  by  0 
atoms  lower  in  the  atmosphere  (refer  to  the  discussion  in  Ref  1  ),  the 
trend  inferred  from  the  data  set  and  one  individual  stable  arc  indicates 
much  more  sensitivity  to  ionization-excitation  altitude  than  is  pre¬ 
dicted  from  a  simple  aerochemical  calculation.  We  recommend  1)  that 
further  individual  flight  data  segments  be  analyzed  to  check  on  this 
important  issue,  and  2)  that  runs  of  more  complete  computer  models  of 
NO  excitation  and  quenching  (0PT1R,  in  particular)  be  made  to  improve 
the  assessment  of  the  expected  (on  the  basis  of  current  DNA  chemical 
reaction  models)  dependenceof  SWIR  yield  on  the  peak  (and  shape)  of  the 
altitude  profile  of  particle  energy  deposition. 

No  statistically  significant  correlation  between  the  SWIR  yield 
and  several  other  auroral  and  measurements  parameters,  in  particular 
the  fractional  solar  illumination,  is  derivable  from  the  segment  of 
data  addressed.  These  parameters  do  not  include  the  amount  of  predosing 
of  the  sampled  region  or  the  temporal -spa t ia 1  stability  of  the  energy 
input;  as  noted  (and  as  further  discussed  in  the  Test  Plan),  the  SWIR 
radiance  lag  would  be  lower  (and  possibly  even  negative)  from  aurora 
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in  which  an  equilibrium  between  energy  input  and  infrared  radiation 
output  is  reached.  The  video  camera  data  contain  the  necessary  informa¬ 
tion  about  the  time-  and  space-variation  of  this  input. 

The  SWIR-emission  output  is  found  to  lag  particle-energy  input, 
with  the  estimates  from  visual  inspection  of  the  radiance  traces  showing 
considerable  variation  among  segments  (6  to  0  sec)  and  no  cases  of 
negative  delay.  We  grouped  the  lag  events  into  two  categories,  >  2  sec 
and  2  sec,  to  investigate  correlation  with  solar  depression  and  frac¬ 
tional  illumination,  auroral  input  intensity  and  peak  height,  type  of 
auroral  input  distribution,  Kp ,  and  measurement  factors  such  as  direction 
of  instrument  pointing  (parallel  to  or  26°  from  the  direction  of  the 
particle  stream)  and  missions  sequence  (potentially  related  to  differ¬ 
ences  in  calibration).  Cast  in  this  form,  the  lag  data  sample  addressed 
do  not  show  sta t i s t i cal ly-s i gn i f icant  dependence  on  any  of  these  para¬ 
meters.  The  effects  of  predosing  and  in  particular  temporal -spat ia l 
stability  of  the  particle  bombardment  (i  .e.  ,  closeness  to  "equilibrium11), 
which  would  tend  to  reduce  the  decorrelation,  have  not  yet  been  con¬ 
sidered;  the  video  camera  data  contain  the  necessary  input  distributions. 

Considering  that  the  ai rcraft-rad iometry  data  represent  the  sole 
source  of  information  on  the  important  issue  of  spat ial -temporal  correla¬ 
tion  of  SWIR  radiation  output  with  energy  input  at  high  (surveillance 
system-scale)  resolution,  it  is  recommended  that  the  complete  data  set 
be  evaluated  and  (in  particular)  that  quantitative  time  delays  be 
derived  from  the  radiance  data  streams.  We  recommend  also  that  high 
altitude  air  chemistry  code  runs  address  specifically  the  question  of 
grow-in  and  decay  of  SWIR  radiance  in  magnetic  zenith-oriented  columns, 
to  validate  the  OoD  SWIR  excitation  model  against  these  simulation  data 
from  the  naturally  disturbed  atmosphere. 
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SECTION  3 


DELAYED  VISIBLE  EMISSION  EXCITED  BY  PRECEDE 


INTRODUCTION 

We  determined  the  yields  of  visible  radiation  after  charged- 
particle  energy  was  deposited  at  E-region  altitudes  by  DNA/AFGL's 
"PRECEDE" rocket  accelerator  (Ref  27  ) .  by  microdens i tomet r ic  analysis 
of  a  photograph  of  the  afterglows.  The  radiance  distribution  resulting 
from  a  series  of  1.0  sec,  0.8  A  pulses  of  2i  kV  electrons  in  the  pre¬ 
vious  ^36  sec  is  shown  in  Figure  13.  The  persisting  self-luminous 
"puffs"  are  centered  in  altitude  at  the  beam-on  positions  along  the 
vehicle's  trajectory  (the  beam  was  off  for  1.0  sec  between  pulses). 

Their  horizontal  offsets  are  due  to  winds,  whose  velocity  perpendicular 
to  the  camera's  line  of  sight  is  as  high  as  200  m/sec  at  some  of  the 
part icl e- i nject ion  altitudes. 

The  photograph  analyzed  was  obtained  by  DNA  contractor  Technology 
International  Corporation  (Bedford,  MA)  who  had  responsibility  for 
processing  the  film  and  relating  its  developed  density  to  log  relative 
exposure.  Its  projection  is  to  the  Tiff  Optical  Site  of  White  Sands 
Missile  Range,  which  is  in  the  vertical  plane  of  the  trajectory  at  a 
mean  slant  range  of  about  115  km.  The  frame  was  exposed  for  A  sec  and 
ends  during  the  73rd  accelerator  pulse,  when  the  rocket  is  at  96  km  alti¬ 
tude  on  downleg  and  moving  vertically  at  just  under  0.7  km/sec.  Camera 
focal  length  was  150  mm  and  relative  aperture  f/0.9  in  the  center  of  the 
field  of  view.  Pulses  73  through  61  are  reasonably  well  defined  in  the 
reproduction,  and  the  "older,"  more  diffuse  glows  from  still  earlier 
energy  input  (pulses  60  through  55,  26  to  36  sec  before  the  exposure  ends) 
are  resolvable  on  the  original  negative. 

DATA  REDUCTION 

Peak  surface  radiances  and  radiances  integrated  over  the  area  of 
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Figure  13.  Photograph  from  the  Tiff  Optical  Site,  WSMR,  of  the  afterglow 

resulting  from  deposition  of  1-sec  pulsed,  2i  keV ,  0.8  A  electrons 
by  PRECEDE.  Exposure  duration  is  4  sec,  f/0.9  on  EK  2485  film. 

It  begins  at  221  sec  after  launch,  when  the  rocket  (bottom  of  the 
bright  streak  containing  pulses  71*73)  is  at  99  km  on  downleg. 

The  vertical  distance  between  centers  of  successive  resolved  "puffs” 
in  the  center  of  the  track  is  about  1  km. 


each  puff  were  determined  by  photographic  photometry,  the  methods  being 
essentially  those  applied  in  Section  1.  The  film  was  scanned  with  a 
75nm-square  microdensitometer  slit,  which  projects  to  65  meters  at  the 
115  km  mean  range  of  the  ser i es  of  el ectron  ejections.  This  is  about  the 
best  spatial  resolution  allowed  by  the  Eastman  Kodak  2485  film's  grain 
noise  and  the  small  increase  in  density  of  the  glow  areas  over  the  sky 
background  density  (which  was  0.26  and  thus  close  to  film  fog  level). 

The  density  increment  was  0.02  in  the  initial  contour  plot  made  by  the 
two-dimensional  scan  of  the  image.  We  removed  stars  and  smoothed  the 
original  film  density  contours  on  the  basis  of  our  experience  with 
similar  calibrated  photographs  of  chemical  releases  in  the  upper 
atmosphere.  The  existing  step  tablet  calibration  was  used  to  assign 
relative  exposures  to  the  contours  in  each  puff. 

Figure  14  plots  the  resulting  maximum  surface  radiances,  the 
minimums  between  these  peaks,  and  the  a rea- i ntegrated  radiances  as  a 
function  of  the  tinv  .nee  the  atmospheric  volume  had  been  initially 
excited  and  partial i,  ionized.  The  data  refer  to  the  4000-6500  A 
sensitivity  range  of  the  camera.  Radiance  integrated  over  puff  area, 
which  we  measured  from  the  isophote  plot  using  a  planimeter,  is  directly 
proportional  to  the  total  rate  of  emission  at  these  wavelengths.  Note 
that  the  maximum  brightnesses  are  comparable  to  the  background  due  to 

O- 

nightglow  (pr inc i pa  1 1 y  ,  the  N0£  continuum),  scattered  starlight, 
zodiacal  light,  and  unresolved  weak  stars;  and  that  the  peaks  are  of 
the  order  of  twice  the  minima  (  in  the  regions  where  puffs  overlap).  The 
maximum  in  peak  radiance  near  10  sec  after  the  initial  disturbance  may 
of  course  be  an  effect  of  the  change  in  energy  deposition  altitude 
rather  than  a  result  of  the  grow-in  and  decay  of  visible  chemilumin¬ 
escence  (the  altitude  range  of  the  measurements  lies  at  and  above  the 
peak  in  the  atmosphere's  oxygen  atom  concentration). 

As  the  density-exposure  transfer  characteristic  of  the  negative 
lacked  an  absolute  radiometric  re'erence,  a  post-calibration  was  needed 
for  quantitative  estimation  of  the  absolute  radiances  and  visible-light 
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yields  from  the  afterglows.  We  applied  both  the  nominal  mean  radiance 
of  the  night  sky  in  the  film's  spectral  sensitivity  range  (Ref  28  )  and 
Eastman  Kodak  data  (Ref  29)  on  its  response  to  derive  an  absolute  ex¬ 
posure  at  the  background  density  of  0.26.  Using  the  night  sky  back¬ 
ground,  an  assumed  lens  transmission  of  0.6  (typical  for  very  fast  lenses 
at  full  aperture),  and  a  Ij-stop  correction  for  2^+85 1  s  reciprocity 
failure  in  the  4  sec  exposure  time  (which  we  had  derived  from  tests  in 

support  of  earlier  chemical  release  imaging  programs),  we  found  that  an 

2 

exposure  of  0.0045  erg/cm  would  produce  this  density.  The  background- 

producing  figure  derived  from  the  film  manufacturer's  data  is 
2 

0.0065  erg/cm  . 

2 

We  adopted  an  absolute-calibration  figure  of  0.0055  erg/cm  at 
the  sky  background,  with  an  estimated  error  of  ±  1 5  stops  or  a  factor 
2.8.  This  error  is  in  large  part  due  to  the  film's  low  contrast  near 
its  threshold  (that  is,  density  is  increasing  very  slowly  with  log 
exposure).  Using  this  normalization  we  calculated  a  total  rate  of  visible 
light  output  of  6(-4,  +11)  watts  from  accelerator  pulse  70  (Figure  14). 
Note  that  the  persisting  yield-rate  from  each  of  the  energy  deposition 
volumes  lies  within  a  factor  2  of  this  figure.  Summing  the  persisting- 
radiation  power,  we  find  that  the  total  yield  in  the  most  aged  puffs 
is  of  the  order  of  100  watt  sec. 

DISCUSSION 

Total  energy  input  to  the  atmosphere  in  each  of  the  accelerator 
pulses  is  2000  watt  sec,  and  the  "prompt"  impact-excited  fluorescence 
efficiency  of  air  within  the  camera's  wavelength  sensitivity  range  is 
less  than  2 X  (assuming  no  neutralization-associated  discharge  takes  place 
in  the  rocket's  vicinity;  none  is  suggested  by  the  groundbased  spectro¬ 
graph  data,  Ref  27).  The  light  yield  during  the  pulse  proper  would 
thus  be  under  40  watt  sec.  It  is  to  say  the  least  surprising  that  a 
comparable  amount  of  energy  (or  more)  comes  out  in  delayed  visible 
radiation;  recall  that  no  similar  afterglow  is  observed  in  the  aurora, 
for  example. 
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No  metastable  species  is  known  to  produce  such  intense  and  long- 
persisting  emission,  and  in  any  case  the  time  dependence  of  the  visible 
light  output  is  inconsistent  with  such  a  simple  decay  process.  The  best- 
known  air  afterglow  mechanism  is  the  reaction  of  nitric  oxide  (known  to 
be  produced  with  high  efficiency  in  irradiated  air)  with  oxygen  atoms, 
which  results  in  the  familiar  yellow-green  glow  (which  is  cyclic  and 
consumes  odd  0) .  A  simple  calculation  shows,  however,  that  to  reach 
the  measured  radiant  intensities  the  concentration  of  NO  molecules  would 
have  to  be  at  least  comparable  to  the  ambient[N2  +  02],  which  is  ener¬ 
getically  not  possible.  No  further  known  process  involving  air  species 
alone  appears  capable  of  producing  this  PRECEDE  afterglow.  It  may  be 
associated  with  effluents  from  the  rocket  motor,  which  were  seen  to 
produce  a  bright  chemiluminescent  wake  along  the  trajectory  (Ref  27); 
if  so,  some  explanation  of  why  the  emission  is  associated  with  the 
ionization-excitation  from  the  electron  pulse  is  needed. 


9^ 


SECTION  4 


EXCEDE:  SPECTRAL  ENERGY-DEPOSITION  DISTRIBUTIONS 


INTRODUCTION 

The  original  function  of  EXCEDE:  Spectral's  ( EX851 .44-1  ,  29  Oct  78; 

Ref  30  )  onboard  radiometric  cameras  was  to  determine  energy  input  rates 
within  the  fields  of  view  of  the  payload's  infrared  spectrometers  and 
radiometers.  The  image  data,  which  are  reviewed  in  Section  2  of  Ref  22, 
are  indeed  directly  suitable  for  this  purpose.  Their  preliminary  analysis, 
however,  shows  that  the  spatial  distributions  of  visible  radiance  differ 
markedly  from  predictions  of  independent  particle-transport  models  (Ref  31) 
both  in  individual  scenes  and  in  their  dependence  on  the  altitude  (i.e., 
ambient  atmospheric  pressure)  at  which  the  3  kV  electrons  are  injected. 

This  finding  is  supported  by  a  recent  analysis  of  some  of  the  groundbased 
camera  data  (Ref  32),  which  disagree  with  the  calculations  at  distances 
from  the  rocket  beyond  the  onboard  camera's  effectively-100  meter  range; 
and  by  the  altitude  profiles  of  air  fluorescence  measured  by  a  photometer 
(Ref  33).  Electric  discharges,  which  are  thought  to  assist  in  neutral¬ 
izing  charge-ejection  spacecraft,  have  a  different  distribution  of 
radiation-exciting  secondary  electrons  than  individual  fast  particles. 

If  the  observed  peculiarity  in  energy  input  pattern  is  due  to  a 
discharge  the  infrared  spectral  intensities  measured  by  EXCEDE's 
instruments  would  differ  from  those  excited  by  energetic  particles  and 
photons  from  atmospheric  nuclear  explosions. 

Hence  emphasis  in  evaluation  of  the  film  data  has  been  directed  to 
determining  volume  emission  rate  distributions,  from  which  glow  excitation 
mechanisms  can  be  inferred.  Feature  intensity  ratios  measured  by  the 
onboard  instruments  also  play  an  important  part  in  determining  the  energy 
spectrum  of  the  electrons  in  the  "plasma"  produced  by  the  ejected  beam 
(into  which  the  instrument  point),  as  is  discussed  in  Ref  22  .  Character¬ 
ization  of  the  ionized-exci ted  volume  near  EXCEDE:  Spectral  applies  not 
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only  in  interpreting  the  measured  infrared  chemi 1 uminous  efficiencies 
but  also  in  designing  spaceborne-accelerator  investigations  in  which  the 
spectrorad iometr ic  measurements  are  made  from  remote  platforms  (the  mother- 
daughter  concept).  The  optical  data  also  provide  estimates  of  the  free 
electron  concentrations,  for  planning  measurements  of  the  infrared  plasma 
(bound-bound)  atomic  radiations. 

Monochrome  film  frames  from  electron-injection  altitudes  1 5~38  km 
below  those  in  Ref  22  have  been  reduced.  The  spatial  and  al t i tude -distri¬ 
bution  data  are  presented  here,  with  a  discussion  of  their  application 
in  inferring  energy  distributions  of  the  secondary  and  degraded-primary 
electrons. 

BACKGROUND 

We  review  briefly  the  photographic  measurements,  particulars  of 
which  are  in  Ref's  22  and  34.  EXCEDE:  Spectral's  30°-conical  electron 
beam  was  directed  upward  along  or  near  the  direction  of  the  geomagnetic 
field,  at  60°  from  the  rocket's  long  axis  and  toward  the  fields  of  view 
of  its  array  of  infrared  instruments  (see  Fig  12  of  Ref  22 ) .  Two  co¬ 
aligned  94°  x  62°-field  cameras,  one  with  color  film  and  the  other  with 
monochrome  film  filtered  for  sensitivity  at  wavelengths  <  4600  A,  viewed 
the  beams  from  a  point  in  or  near  the  meridian  plane  3i  meters  south  of 
the  accelerator  segment  that  provided  most  of  the  useful  data.  In  this 
projection,  the  field  line  onto  which  the  electrons  are  injected  lies  inside 
the  camera's  image  field  from  "infinity"  down  to  within  1i  meters  from 
the  port  from  which  the  electrons  emerge.  The  cameras'  viewing  perspective 
provides  good  lateral  spatial  resolution  to  at  least  30  meters  from  the 
rocket  (and  useful  data  to  MOO  m)  ,  while  the  first  refocusing  node  of 
3  keV  electrons  in  the  geomagnetic  field  over  Poker  Range  (where  the 
rocket  was  launched)  is  at  21  meters.  No  node  appears  to  be  above  detec¬ 
tion  threshold  in  any  of  the  500-odd  data  frames  (although  one  may  be 
present  in  an  early  frame  in  which  <  0.1  amperes  was  ejected). 

Sequences  of  8,  80,  and  800  millisec  exposures  were  made  to  extend 
the  dynamic  range  of  scene  brightness  that  could  be  recorded.  The 
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black-and-white  film  is  sensitive  to  N*  First  Negative  and  Second 
Positive  band  air  fluorescence,  and  insensitive  to  the  persisting  5577  A 
01  line.  The  onboard  cameras'  high  resolution  close-in  to  the  rocket  — 
in  particular,  in  the  region  into  which  the  infrared  instruments  point  -- 
is  complemented  by  lower-resolution  (>  30  m  )  photographs  of  the  complete 
beam  energy  deposition  pattern  from  groundbased  cameras. 

The  films,  which  are  the  first  from  onboard  a  charge  ejection 
vehicle,  show  a  field  aligned  "core"  with  the  magenta  halo  most  probably 
due  to  excitation  of  N2  triplet  state  transitions  (First  and  Second 
Positive  bands)  by  slower,  laterally-scattered  secondary  electrons.  Several 
of  the  radiance  distributions  on  the  monochrome  negative  have  been  re¬ 
duced  manually,  starting  from  equi-density  contour  plots  made  with  a 
two-dimensional  1 y- scanning  microsens i tometer  (as  in  Sections  1  and  3). 

The  brightness  distributions  when  7  amperes  was  ejected  at  97.7  and 
8^.9  km  on  downleg  are  shown  in  Fig  15(lens  vignetting  has  been  corrected 
out).  Similar  plots  for  123  km  upleg  and  113  km  downleg,  with  further 
clarification  of  the  camera  perspective  and  a  plot  indicating  the  hori¬ 
zontal  scale,  are  in  Ref  22.  These  "maps"  serve  as  starting  points  for 
tracing  radiance  along  phys ica 1 1 y-descr i pt i ve  directions  such  as  the 
core  axis  (Fig  16) ,  and  for  unfolding  to  generate  volume  emission  rates  -- 
the  "fundamental"  indicator  of  plasma  conditions  present  in  the  images. 

DATA 

As  Fig's  15  and  16  show,  within  a  few  rocket  diameters  the  glow's 

surface  brightness  increases  very  rapidly  toward  the  accelerator  and 

rocket  body.  Beyond  this  range  the  increase  with  distance  is  due 

principally  to  the  increasing  sight  path  through  the  (optical ly-thin) 

beam,  that  is,  to  viewing  perspective.  (The  minimum  is  also  present  in 

Fig's  20  and  22  of  Ref22.)  The  glow  has  about  twice  the  initial  angular 

divergence  measured  in  laboratory  tests  with  the  electron  accelerator 
-4  -3 

operated  into  10  -  10  torr  air  pressure,  and  as  noted  no  effects  of 

refocusing  show  in  the  five  isophote  plots  (the  one  for  116  km  upleg, 
not  reproduced  here,  qualitatively  resembles  the  others). 
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The  radiance  distribution  along  the  beam's  field-aligned  axis, 
both  in  the  region  close-in  —  whose  "shape"  is  strongly  at  variance  with 
the  independent-particle  transport  model  --  and  beyond  the  minumum  changes 
little  (and  not  even  monotonica I ly  with  altitude)  over  the  factor-250 
range  of  ambient  air  density.  While  the  close-in  falloff  might  be  ex¬ 
plained  by  invoking  a  local  discharge,  or  rocket  outgassing  (which  is 
known  to  play  a  part  in  discharges,  Ref  35),  the  small  (factor  <  3) 
variation  with  altitude  at  10  -  100  meters  from  the  rocket  would  not  appear 
to  be  explained  by  an  outgassing-related  process. 

A  narrow  range  of  emission  intensities  was  also  measured  by  an 
onboard  photometer  pointed  near  the  minimum-radiance  area  on  the  films 
(Ref  33).  This  instrument  showed  much  high  signals  on  upleg  than  on 
downleg  at  the  same  altitude.  The  effective  ^-molecule  concentrations, 
assuming  that  this  increase  is  due  to  outgassing  (only),  are  given  at 
top  right  of  rig  16.  Additionally,  ground  camera  data  show  that  the 
energy  deposited  at  123  km  ejection  altitude  does  not  fit  an  independent- 
particle  model  in  either  the  first  few  hundred  meters  from  tha  vehicle 
(where  it  is  too  high)  or  in  the  )  to  3  km  region  (too  low)  (Ref  32). 

This  finding  indicates  a  redistribution  of  the  energy  of  the  initial 
injected  beam,  beyond  that  resulting  from  Coulomb  collisions  with  atoms 
and  molecules  of  ambient  atmosphere.  We  note  that  the  ground  photo¬ 
graphy  (at  the  one  ejection  altitude  analyzed,  123  km)  shows  excess 
emission  to  ranges  well  beyond  the  ^5  meter  minimum  in  Fig  16  . 

DISCUSSION 

Models  of  infrared  surveillance-systems  performance  implicitly 
assume  that  the  particles  output  from  nuclear  explosions  --  thermal 
x-rays,  debris  ions/atoms,  6  and  y-rays  —  act  independently  of  one 
another  in  exciting  the  distant  atmosphere  (except  when  the  doses  are 
high  enough  to  produce  substantial  heating).  If  the  distribution  of 

energies  of  the  secondary,  tertiary,  .  electrons  in  the  air  volume 

probed  by  EXCEDE  differs  from  that  from  nuclear  particles,  the  relative 
populations  of  infrared-emitting  species  (and  their  precursors)  will 


10  Phot°graphic  Column  Brightness/  Erg/cm  -sec-sr 


Figure  16.  3800-4600  A  EXCEOE:  Spectral  radiance  along  the 

injected  electron  beam's  axis,  projected  to  the 
onboard  camera.  Ambient  Nj  density  at  the  five 
ejection  altitudes  spans  a  range  of  about  250. 
Also  shown  are  the  radiances  measured  by  a 
3914  A-band  photometer  pointed  near  the  glow 
minimum. 
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be  different,  which  would  impair  the  accuracy  of  EXCEDE's  simulation. 

For  example  an  excess  of  electrons  of  kinetic  energy  near  3  eV  would 
result  in  enhanced  emi ss ions  from  CO^  molecules,  both  by  direct  impact 
excitation  and  (in  particular)  increased  energy  transfer  from  mole¬ 
cules,  whose  cross-section  for  vibrational  excitation  peaks  near  this 
electron  energy. 

This  effect  of  this  energy  distribution  on  some  of  the  atmosphere's 
emissions  is  shown  schematically  in  Fig  17.  After  a  few  ionization  mean 
free  paths,  an  i ni t ial 1 y-monoenerget ic  electron  beam  is  accompanied  by  a 
comparable  number  flux  of  slower  secondary  electrons  and,  because  of  the 
statistics  of  collisions,  becomes  spread  in  energy.  (The  energy  spectrum 
shown  in  the  top  panel  is  in  effect  the  DPOSIT  distribution,  which  is 
input  to  ROSCOE's  model  for  air  fluorescence  (Ref  36).)  These  secon¬ 
daries  vibrational  ly  excite  N2  and  electronically  excite  N^s  triplet 
manifold,  which  --  as  pointed  out  in  Ref  22  --  was  unexpectedly  strong 
in  EXCEDE's  spectrums.  Additionally,  non-homopolar  molecules  such  as 
CO^  (and  H^O ,  whose  infrared  radiations  were  also  strong  at  EXCEDE)  have 
high  cross-sections  for  vibrational  excitation  by  soft  electrons. 
Discharges  produce  large  numbers  of  secondaries  with  energies  of  a  few 
eV  (in  comparison  to  impact- i on i zat ion  by  the  primary  electrons  (see 
Fig  17a),  as  shown  in  both  laboratory  and  rocket  experiments  (Ref  37). 

AFGL/OPR  is  currently  evaluating  spectroscopic  data  from  EXCEDE: 
Spectral  (Ref  33),  to  determine  their  fit  to  laboratory  and  auroral- 
measurement  values.  (Auroral  electrons  effectively  simulate  nuclear  - 
particle  excitation,  as  they  have  essentially  the  same  secondary  electron 
spectrum.)  The  ratio  of  column  intensities  in  N2  Second  Positive  to  N+ 
First  Negative  bands  is  indicated  in  Fig  1 8 ,  which  shows  quantitatively 
the  effect  of  MO  eV  electrons  on  N2  triplet  and  vibrat  ional -state 
exci tat  ion. 

Discharges  of  two  principal  types,  static  (in  effect  corona,  Ref  38) 
and  beam-plasma  (Ref  39) ,  have  been  invoked  as  the  source  of  the  return 
current  needed  to  neutralize  the  electron-ejection  vehicle  and  beam.  The 
flux  of  ambient  electrons  from  the  natural  ionospheric  plasma  is  only  of 
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Figure  17.  Model  secondary  electron  flux  spectra  due  to  independent  primary- 
particle  collisions  and  to  discharges  (a),  and  cross-sections  for 
impact  excitation  of  N  by  1-1 04  eV  electrons  (b)  . 
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Figure  18.  Excitation  cross-sections  and  energy  spectra  of  thermal 
electrons  at  temperature  0.87  -  17.5  eV.  The  spectrum 
measured  on  charge-ejection  rocket  Echo  III  (Ref  37) 
fits  no  unique  temperature.  Ratios  of  emission  intensity 
in  the  features  indicated  are  in  parentheses. 
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the  order  of  1  milliamp  per  m  of  vehicle  surface  (Ref  39  )  ,  and  the 
departing  electron  beam  is  not  believed  to  produce  enough  new  free 
electrons  nearby  to  maintain  its  potential  at  the  observed  (Ref  38  ) 

10's  of  volts.  Some  observables  of  these  discharges  are  listed  in 
Table  8.  The  static  discharge  is  driven  by  the  space  plasma-shielded 
vehicle  potential  ,  which  produces  a  return  current  of  MOO  eV  electrons 
(which  is  near  the  maximum  of  the  excitation  cross-section)  that  create 
still  further  secondaries  near  the  vehicle.  Beam-plasma  interactions 
result  in  an  RF  discharge,  in  which  can  exist  electric  fields  that  broaden 
the  energy  distribution  of  the  initial  beam  electrons  themselves  (as 
indicated  in  Fig  17a)  and  so,  among  other  effects,  destroy  the  refocusing. 

It  should  be  noted  that  simple  Coulomb  self-repulsion  of  the  7 
ampere  ejected  beam  could  also  be  causing  its  broadening  and  apparent 
failure  to  refocus,  as  its  expected  "egui 1 ibrium"  radius  in  a  partially- 
ionized  medium  can  be  several  meters  (Ref  39).  This  space-charge  effect 
could  smear  out  the  node  so  that  its  radiance  contrast  falls  below  de¬ 
tectability  in  the  photographic  projections. 

The  spatial  distribution  of  volume  rate  of  excitation  can  be 
derived  by  unfolding  radiance  distributions  such  as  are  shown  in  Fig  15, 
as  the  glow  appears  to  a  good  approximation  cyl indrical ly  symmetric. 

The  ionization-dissociation  rate  is  proportional  to  the  rate  of  excita¬ 
tion  of  3911*  A  and  A278  A  N+  First  Negative  band  radiation,  which  is  the 
principal  component  of  the  photographic  signal  (the  Nj  Second  Positive 
bands'  contribution  can  be  estimated  from  the  spectrum  data).  Summing 
of  these  volume  emission  rates  thus  determines  the  total  energy  input 
rate  within  the  fields  of  individual  onboard  instruments. 

These  exc i tat  ion- ion i zat ion  rates  also  serve  in  estimating  free 

electron  densities,  which  are  needed  to  assess  the  intensities  of  plasma 

emissions.  At  equilibrium  between  creation  and  chemical  loss  of  electrons, 

and  for  electron  densities  n^  in  the  expected  range  (where  two-body 

reactions  dominate),  Q  =  an  ^  ,  where  a  is  the  recombination  coefficient. 

e 

Dissociative  recombination  with  molecular  ions  is  the  dominant  neutral iza- 
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Table  8 . 

Basic  Concept 

Current/Voltage 
at  M40  km  al  t  i  tude 

Threshold 

Effect  on  initial  beam 
radiance  distribution 


Opt ica I  f I i cker 

RF  emission 

Secondary  electrons 


Characteristics  of  "corona"  discharge  and 
Beam  Plasma  discharge 


Corona 


BPD 


DC  avalanche  of 
ambient  electrons 
driven  by  electro¬ 
static  field  from 
vehicle 


Instability  in  ambient 
plasma  driven  by  beam's 
charge,  RF  fields  grow 
to  accelerate  electrons 
above  the  ionization 
threshold  of  atmospheric 
molecules /atoms 


M  mA,  3  kV 


%80  mA,  3  kV  threshold 


Continuous  process  Sharp  but  with  some 

structure  near  onset 


Vi rtuel ly  none  i f 
the  return  current 
is  isotropic;  within 
a  few  rocket  radii 
if  the  current  flows 
in  the  path  ionized 
by  the  ejected  beam 


Pitch  angle  scattering  re¬ 
sults  in  lateral  broadening, 
intensified  backward  beam 
component,  focusing  node 
destruction;  factor-10  in¬ 
crease  in  total  fluorescence 
yield  within  ^200  meters; 
longitudinal  distribution 
warped  by  heating  of  beam 


Not  observed,  may 
be  present 


At  ion  cyclotron  frequency 
^20  Hz 


No 


At  and  above  plasma  fre¬ 
quency,  complex  spectrum 


It  is  expected  that 
distribution  can  be 
probably  results  in 
at  each  energy,  and 
of  ^100  eV  electror 
primaries. 


some  differences  in  the  energy 
identified  in  the  theory.  BPD 
a  broader  angular  distribution 
should  show  "suprathermal  tails" 
which  behave  much  like 
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tion  reaction,  when  the  density  of  molecular  ions  is  >  10  ^  that  of 
atomic  ions,  as  would  be  expected  under  EXCEDE:  Spectral's  conditions. 
(This  is  the  reverse  of  using  free  electron  densities  measured  by  the 
Chatanika  incoherent-scatter  radar  for  determining  energy  input  rates  in 
the  auroral  ionosphere,  which  also  requires  knowledge  of  a; see  Ref  40.) 

We  note  that  estimates  of  n  in  discharges  produced  by  electron  injection 

8  e  9  3 

vary  from  less  than  10  to  more  than  10  per  cm  . 

The  spatial  distributions  of  volume  emission  rate  that  result  from 
the  neutralization-assisting  processes  would  have  the  following  general 
propert i es . 

Corona:  varying  relatively  slowly  with  distance 
from  the  rocket  (Ref  38),  and  not  confined  to 
the  initial  beam  volume. 

Beam-plasma  discharge:  lateral  broadening  due 
to  pitch  angle  scattering  of  the  ejected  electrons, 
and  an  axial  distribution  altered  by  heating  of  the 
beam  (that  i s  ,  i t  1  s  spread  to  both  higher  and  lower 
energies);  confined  to  the  (expanded)  beam  volume. 

We  note  that  beam-plasma  discharge  is  not  expected 
to  occur  at  the  air  pressures  below  ^95  km  altitude; 
and  that  the  film  data  so  far  do  not  indicate  any 
sharp  transitions  in  the  radiance  pattern. 

Outgassing:  confined  to  the  region  partially 
ionized  by  the  outgoing  beam;  depends  on  the 
distribution  in  space  of  the  outgassed  molecules, 
which  can  be  calculated. 

As  noted,  identification  of  the  neutralization  mechanism  would  provide  an 
indication  of  the  effect  of  secondary  electrons  on  EXCEDE:  Spectral's 
infrared  data.  The  primary  effect  of  outgassing  is  to  change  the  apparent 
altitude  profiles  of  the  infrared  yields;  while  the  emission  spectrum 
may  be  changed  by  the  introduction  of  new  chemical  species  (H^O  in 
particular),  outgassed  molecules  do  not  change  the  flux  of  secondary 
electrons  relative  to  the  primary  beam. 
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SECTION  5 


RADIANCE  STRUCTURE  OF  THE  LOWER  THERMOSPHERE 


BACKGROUND 

Air  Force  Geophysics  Laboratory  recently  made  a  six-month  satellite 
survey  of  the  spectral  intensity  of  1100  -  2900  A  ultraviolet  radiation 
from  the  earth's  atmosphere,  with  high-resolution  spectrometers  and  a 
wavelength-cycling  narrow-field  photometer  on  a  polar-orbiting  vehicle 
(Ref  41).  The  spacecraft,  DoD  S3_A,  was  near  sun-synchronous  at  an 
inclination  of  about  96.5°  and  crossed  the  equator  at  1030  and  2230  LT. 

The  radiance-measuring  instruments  (of  experiment  CRL-726)  viewed  directly 
toward  the  nadir  from  altitudes  between  160  and  260  km. 

When  the  thermosphere  is  sunlit,  almost  all  of  its  vacuum  UV 
emission  (wavelength^  1800  A)  is  fluorescence  of  N^  and  0  excited  by 
photoelectrons  ejected  from  air  molecules  and  atoms  by  extreme  UV  solar 
photons  (wavelength  .$  600  A).  Absorption  by  the  atmosphere's  0^  molecules 
in  their  Schumann-Runge  continuum  prevents  most  VUV  photons  that  originate 
below  approximately  110  km  (  the  exact  number  depends  on  the  wavelength) 
from  reaching  satellite  altitudes.  Variations  in  the  mixing  ratios 

would  vary  the  radiances  of  the  N^  and  0  dayglow  features, 
through  the  change  in  both  the  concentrations  of  target  (emitting)  species 
and  the  energy  spectrum  of  the  secondary  electrons, and  also  because  of 
the  change  in  attenuation  of  upward-directed  VUV  photons.  It  is  safe 
to  assume  that  the  incident  EUV  flux's  variation  has  only  small  fre¬ 
quency  components  at  the  spatial  scales  of  surveillance-systems 
interest  (a  few  km,  which  considering  the  satellite's  orbital  velocity 
of  8  km/sec  translates  to  MHz  in  the  solar  irradiance).  Thus  fluctua¬ 
tions  in  the  VUV  column  intensities  measured  from  S3_i»  would  provide  a 
measure  of  the  atmospheric-species  structure  at  altitudes  from  which 
this  radiation  is  emitted. 
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Spatial  structure  of  typical  scale  size  30  km  near  150  km  altitude 
has  been  observed  at  high  latitudes  by  satellite-borne  probes  (Ref's 
42  ,43  ).  The  variations  in  the  concentration  of  argon  atoms  and  in 
mixing  ratios  are  associated  with  geomagnetic  storms.  These  in  s i tu 
measurements  suggest  a  similar  inhomogeneity  of  infrared-active  molecular 
species  such  as  CO^,  which  would  lead  to  irregularities  in  the  sky  back¬ 
ground  presented  to  IR  surveillance  sensors  when  the  upper  atmosphere 
absorbs  energetic  radiations  from  nuclear  explosions;  this  "clutter" 
mimics  radiating  targets  and  otherwise  degrades  system  performance. 

The  VUV  radiance  distribution  remotely  sensed  from  S3-4 
provides  an  opportunity  to  investigate  further  this  variability  of 
species  concentrations  in  the  lower  thermosphere,  particularly  at  high 
spatial  frequencies  (>,  1  cycle/km,  the  photometer's  footprint  being 
<  1  km).  Hence  we  have  made  a  preliminary  analysis  of  the  photometer 
data  from  one  flight  segment,  to  determine  if  the  nadir  radiance  distri¬ 
bution  differs  from  that  expected  and,  if  so,  to  assess  the  effort 
required  to  quantify  the  spatial  inhomogeneity. 

DATA  BASE 

S3”4's  photometer  isolated  M50  A  FWHM  wavelength  bands  centered 
near  1216,  1340,  1550,  and  1/ 50  A,  with  full  field  of  view  5.9,  1.65, 

0.46,  or  0.12  degrees  circular  (details  are  in  Ref  41  ) .  An  examination 
of  published  (in  Ref  41  )  daytime  data  traces  at  those  altitudes  where 
the  radiance  signals  are  virtually  uncontaminated  by  aurora  indicated 
that  the  1550  A  data  have  the  largest  short-period  (MO  sec,  80  km 
satellite  traverse)  fluctuations.  This  finding  would  be  expected 
because  the  photometer  channel  measures  principally  the  optical ly-thin 
a1^-  X  ^g+  Lyman-Bi  rge-Hopf  ield  bands  (largely  in  the  Av  =  - 1  to -4 
sequences),  which  are  emitted  from  altitudes  135  to  180  km  FWHM  (Ref  44  ) . 
Smooth-appearing  data  traces  are  shown  from  the  1340  A  channel  (which  responds 
principally  to  the  optical ly-thick  01  A  resonance  triplet),  the  1216  A 
channel  (mostly  opt ical ly-thick  Lyman-a  along  with  some  01  resonance 
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radiation),  and  in  the  1750  A  channel  (whose  signal  is  dominated  by 
Rayleigh  scatter  from  altitudes  well  below  100  km). 


We  therefore  selected  for  this  first  analysis  a  segment  of  the 
daytime  radiance  from  this  photometer,  whose  half-maximum  response  is 
at  1485  -  1650  A  and  tenth  -  maximum  at  1470  -  1740  A.  Its  field  of 
view  had  been  set  at  1.65°  during  the  segment  chosen.  The  data  from 
the  complete  satellite  revolution  (445,  13  Apr  78),  compressed  and  pro¬ 
cessed  to  nadir  radiance  (and  hand-smoothed),  are  in  Figure  7  of  Ref  41 . 
Figure  19  is  a  plot  of  the  number  of  counts  in  1/100  of  the  0.0100-sec 
sampling  intervals,  in  the  period  when  the  solar  zenith  angle  is  less 
than  100°  (so  that  the  thermosphere  is  sunlit).  The  full  set  of  data 
points  for  this  period  was  provided  to  PhotoMetrics  by  AFGL  as  a  numer¬ 
ical  listing  and  on  magnetic  tape  in  a  format  suitable  for  processing 
by  AFGL's  CDC-6600  computer. 

Figure  19  g ives  in  addition  to  universal  time  the  local  (solar)  time 
and  solar  zenith  angle  at  the  VUV  instrument's  lat i tude- long  i  tude  ,  which 
are  listed  along  with  its  altitude.  We  narrowed  the  initial  investigation 
down  to  a  sunlit  region  at  high  latitude  but  well  south  of  the  auroral 
oval.  Specifically,  we  chose  the  ^23,000  VUV  radiance  data  points 
accumulated  between  0428:01  and  0431:52.9,  when  the  satellite  was  moving 
southward  between  77i°  and  62°N  latitude  (over  north  central  Siberia). 
During  these  four  minutes  the  instrument's  altitude  decreases  from 
17H  to  166  km,  that  is,  it  is  moving  downward  throug.  the  L-B-H  band- 
emitting  altitudes.  As  its  latitude  decreases  the  sun's  elevation 
angle  --  which  determines  the  EUV  flux  incident  on  the  atmosphere  -- 
increases,  from  21°  to  36° •  The  mean  VUV  radiance  is  seen  to  increase 
by  a  factor  1.63  during  the  data  segment,  while  the  component  of  solar 
irradiance  normal  to  the  horizontal  increases  by  a  factor  1.64  (=  sin  36V 
sin  21°).  Hence  we  may  conclude  that  the  increase  in  EUV  flux  has  a 
much  greater  effect  on  the  atmosphere's  radiance  than  the  decrease  in 
absorption  of  this  energy  in  the  column  defined  by  the  nadir-pointing 
instrument's  field.  The  ^5  sec-averaged  radiance  increases  about  13% 
per  min  and  the  "footprint"  of  the  photometer's  field  of  view  at  the 
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Figure  19-  Number  of  1550  A  photometer  counts 
interval  in  the  daytime  section  of 
Satellite  S3-1*.  The  data  segments 
dicated  at  upper  right. 
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altitude  of  peak  emission,  150  km,  decreases  from  620  to  460  meters 
over  the  data  period;  the  spacecraft's  8  km/sec  orbital  velocity  moves 
this  field  horizontally  one  projection  diameter  in  0.08  to  0.06  sec. 

Plots  in  Ref  45  of  the  Q  =  3  auroral  oval  indicate  that  the  segment 
of  satellite  track  is  >  1000  km  south  of  the  oval;  Kp  0=Q)  was  only  2+  at 
the  time  of  measurement  (Ref  46).  In  fact,  the  northern  oval,  being  centered 
near  280°  longitude,  77°  latitude,  rarely  expands  below  77i°N  at  the  satel¬ 
lite's  110°  longitude.  The  absence  of  auroral  contamination  of  the  VUV 
radiance  is  further  qualitatively  indicated  by  lack  of  obvious  spikes 
in  the  trace.  The  magnetic  activity  index  Ap  for  11,12,  and  13  Apr  78 
was  64,  29,  and  24  (Ref  46);  the  most  magnetically  disturbed  day  of  the 
month  was  two  days  before  the  13  April  data  period. 

DATA  ANALYSIS 

In  the  absence  of  geophysical  or  instrumental  perturbations,  the 
number  of  counts  of  random  independent  events  in  sampling  intervals  of 
fixed  duration  will  be  Poisson-distributed.  Since  the  average  number  of 
1550  A-band  photons  registered  each  1/100  sec  is  as  high  as  100,  the 
Poisson  distribution  can  be  extremely  closely  approximated  by  a  normal 

(Gaussian)  whose  variance  is  equal  to  the  mean  (see,  for  example  Ref  47). 

-1  _  z  _ 

That  is,  the  probability  density  at  N  counts  would  be  (2tt)  2  exp -(N-N) /2N , 
where  N  is  the  mean  or  expectation  number  (about  100,  as  noted). 

Instrument  instabilities  —  in  photometric  response,  or  as  jitter 
in  the  length  of  the  sampling  interval  —  were  found  by  the  experiment 
group  to  have  a  negligible  effect  on  the  seconds-sca le  variations  in 
count  rate.  Dark  current  averages  less  than  1/100  count  per  sampling 
interval,  and  its  noise  in  any  case  adds  in  quadrature  to  that  due  to 
the  natural  fluctuation  in  rate  of  arrival  of  VUV  photons. 

The  geophysical  effects  on  the  atmosphere's  VUV  radiance  include 
the  long  period  change  due  to  the  increasing  elevation  angle  of  the  sun, 
evident  in  Figure  19;  variations  due  to  variable  EUV  output  from  the  sun, 
which  as  noted  are  small  at  the  time  scales  of  interest;  and  the  afore- 
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mentioned  i nhomogene i ty  of  species  concentrations  limits.  Per  iodic i t ies  in  the 
radiance  data  stream,  which  would  appear  as  peaks  in  its  power  spectrum, 
would  help  identify  the  source(s)  of  the  atmospheric  perturbation.  For 
example,  spatial  structure  might  be  found  to  have  characteristic  wave¬ 
lengths  comparable  to  those  of  traveling  ionospheric  disturbances. 

Considering  that  the  magnitude  (or  even  existence)  of  fluctuations 
beyond  those  expected  from  pure  photoelectron  counting  statistics  is 
not  established,  however,  we  judged  it  more  cost-effective  to  first 
apply  simple  statistical  tests  (rather  than  Fourier-transform  methods) 
to  the  data  set. 

We  applied  (from  the  large  number  of  available  tests  of  the  statis¬ 
tics  of  data  sets)  Student's  t-test  to  determine  whether  the  mean 
and  variance  are  equal  within  statistically-allowed  limits  and  the  chi- 
square  test  to  determine  how  well  the  daytime  VUV  data  fit  the  afore¬ 
mentioned  normal  distribution.  It  was  necessary  to  filter  out  the  effect 
of  the  slow  trend  in  tne  signal  ,  which  in  fact  turned  out  to  reduce  the 
reliability  of  the  analysis. 


Student's  t~test 

This  well-known  procedure  tests  a  hypothesis  about  the  mean  of  a 
sample  of  data.  We  hypothesize  first  that  the  photoelectron  counts  are 
Poisson  distributed,  that  is,  that  their  mean  is  equal  to  the  variance. 
The  parameter  t  is  defined  as 


_ (mean  of  sample)  —  ("true"  mean) _ 

(standard  deviation  of  sample)/(number  of  data  points 


n) 


i  ’ 


The  probability  that  t  is  less  than  p,  or  percentile  of  the  t  distribution 
for  the  number  of  degrees  of  freedom  (which  is  to  an  excellent  approxi¬ 
mation  n) ,  is  2p  -  1.  That  is,  when  t  <  p  we  are  confident  at  the  2p  -  1 
level  that  the  population  (true)  mean  lies  within  (standard  deviation)  x 
( t//n)  of  the  sample  mean.  If  the  standard  deviation  of  the  VUV  radiances 
lies  between  these  confidence  limits  we  accept  the  hypothesis  that  they  are 
Poisson  distributed,  which  is  to  say,  that  the  photons  are  arriving  at  the 
satellite  as  independent  quantum  events  from  a  uniformly-emitting  atmosphere. 
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As  the  low-frequency  trend  in  the  data  has  the  effect  of  increasing 
the  variation,  we  attempted  to  remove  the  trend  by  filtering.  We  cal¬ 
culated  running  means  of  length  £  consecutive  samples  --  low  pass 
filtering  --,  which  when  subtracted  from  each  individual  data  point  -- 
high-pass  filtering  --  makes  the  new  mean  very  close  to  zero  (detrending). 
Table  9  shows  results  of  application  of  these  filters  with  £  up  to  300 
on  segments  of  the  data  stream. 

Without  detrending  (£  =  1)  the  variances  for  sample  lengths  greater 
than  500  data  points  (5  sec  of  satellite  data)  become  unrealistically 
high,  and  with  £  =  20  (l/5-sec  integration)  they  remain  considerably 
larger  than  the  means.  Application  of  the  t-test  with  95%  and  99*9% 
confidence  limits  (Table  lo)  shows  the  variance  tobejust  outside  the 
expected  range  when  the  complete  2300-point  data  set  is  high-pass-filtered 
with  £  =  50;  it  also  outside  the  range  in  two  of  the  three  subsets  with 
£  =  20.  This  means  that  the  hypothesis  of  a  Po i sson-d i st r  i bu ted  set  of 
VUV  radiance  data  is  rejected  at  these  levels  of  confidence. 

Chi-square  Test 

This  test  (Ref  i»8 )  assesses  the  "goodness  of  fit"  of  experimental 

data  to  a  preselected  statistical  distribution.  In  light  of  the  failure 

of  the  assumption  of  a  Poisson  distribution,  we  adopted  the  hypothesis 

that  the  radiance  readings  are  Gaussian  (normally)  distributed.  To 

minimize  the  effect  of  the  long-period  trend,  we  divided  the  data  set 

into  3-2  sec  intervals,  during  which  the  mean  count  rate  changes  by 

<  i%;  and  we  treated  the  first  20  such  contiguous  internals  after 

0A28:01.  The  procedure  involves  finding  the  mean  and  standard  deviation 

of  each  data  subset,  placing  the  readings  into  K  bins  of  uniform  count- 

number  width,  and  then  calculating  *  £  (f  -  f  )^/f  ,  where  f  is  the 

I  O  P  6  O 

observed  occurrence  of  the  number  of  counts  and  f  the  expectation 
occurrence  or  frequency  (in  terms  of  a  normal  distribution  having  the 
computed  mean  and  standard  deviation).  We  applied  standard  procedures 
for  lumping  together  the  sparse  number  of  data  points  that  lie  several 
standard  deviations  from  the  mean,  which  improves  the  statistics. 
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lacle  3,  High  and  low- pass  filtering  of  the  S3 
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Table  10.  Results  of  t-test  of  the  S3 -4  VUV  data*. 


Number  of 
Sample 
Points 

Filter 

Length 

i 

Mean 

N 

Variance 

95% 

Confidence 

Limits 

99.  9% 
Confidence 
Limits 

500 

20 

83.39 

79.61 

±0.782 

± 1. 313 

1150 

20 

86.75 

87. 18 

±0.  539 

±0.906 

2300 

20 

88.73 

85.67 

±0.  378 

±0.635 

23295 

50 

115.91 

117.23 

±0.  139 

±0.233 

95% 

Confidence 
Limits  on  Mean 

Variance  i 
Within  j 
Limits  ?  i 

99.9% 
Confidence 
Limits  on  Mean 

Variance 
Within 
Limits  ? 

500 

82.6  -  84.17 

1 

No  ! 

| 

82.07 

-  84.7 

No 

1150 

86.21  -  87.28 

Yes  ! 

| 

85.84 

-  87.65 

Y  es 

2300 

88.35  -  89.10 

No  ! 

| 

88.09 

-  89.36 

No 

23295 

115.77  -  116.04 

No  ! 

115.67 

-  116.14 

No 

$ 

High  pass  filtered;  starting  042801,  revolution  445 
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The  number  of  "degrees  of  freedom"  for  the  test  is  taken  as  K-2 

because  the  two  parameters  mean  and  standard  deviation  are  taken  from 
2 

the  data  sample,  x  is  tabulated  (in  Ref's  47  and  48  ,  for  example)  as 

a  function  of  the  number  of  degrees  of  freedom  and  the  statistical 

2 

significance  level.  If  the  x  calculated  from  the  experiment  data  set 
is  less  than  the  tabulated  value,  the  hypothesis  is  acceptable  at  the 
significance  level  selected. 

We  adopted  for  computational  convenience  a  bin  size  of  b  photo¬ 
electron  counts  (for  example,  96  to  99  in  the  0.01  sec  interval)  after 
preliminary  calculations  indicated  that  the  result  was  not  critically 
dependent  on  this  figure.  The  resulting  means,  variances  (square  of 
standard  deviation),  and  the  chi-squares  for  the  99-5%  significance 
level  are  shown  in  Table  11.  It  is  found  that  19  of  the  20  data  subsets 
(all  except  the  10th)  are  in  agreement  with  the  hypothesis  of  a  normal 
distribution  at  this  significance  level.  This  finding  provides  a  very 
strong  indication  that  the  VUV  radiance  data  in  this  6*4  sec  measurement 
interval  are  normally  distributed. 

The  variances,  however,  are  somewhat  larger  than  expected  from 
pure  photon-arrival  statistics,  which  suggests  that  some  further  random 
process  is  taking  place.  13  of  the  20  variances  are  greater  than  their 
corresponding  means  and  the  mean  of  the  variances  (Av  1-20  entry  in  Table 
11  )  is  bi%  greater  than  the  mean  of  the  means.  This  finding  implies 
that  a  phenomenon  other  than  photon  noise  could  be  impacting  the  radi¬ 
ance  data.  Its  rms  magnitude  would  be  about  30%  of  the  M0%  (100  ^) 

2  i 

relative  photon  noise  fluctuation  (as  1.0^5  =  (l  +  0.3  )  ),  or  3%  of  the 
mean  signal  amplitude. 

SUMMARY 

As  noted,  the  statistical  tests  applied  to  CRL-726's  VUV  data  are  in¬ 
tended  to  determine  whether  fluctuations  in  the  nadir  radiance  measured 
from  the  S3-**  satellite  over  and  above  quantum  (photon,  or  more  strictly 
speaking  photoelectron)  noise  are  present.  A  strongly  positive  result 
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Table  11. 

Ch i -square 

Test  of  VUV 

Data 

Data 

Mean 

Sect  ion 

K 

Counts 

V a r i ance 

2 

X 

1 

8 

86.8 

83.I 

16.44 

2 

10 

88.1 

104.0 

10.44 

3 

9 

89.2 

92.6 

18.70 

4 

9 

89.6 

99.2 

9.98 

5 

9 

89.3 

82.8 

21 .64 

Tabulated 

6 

9 

90.9 

93.5 

17.44 

2 

X 

7 

9 

92.0 

84.3 

17.05 

for 

99.5% 

8 

10 

93-3 

111.8 

14.68 

K 

2 

X 

9 

10 

92.7 

104.5 

14.34 

8 

22.0 

10 

10 

93.9 

112.6 

28.12 

9 

23.6 

1 1 

9 

3b.  0 

96.3 

14.23 

10 

25.2 

12 

10 

96.0 

91.9 

22.26 

1 1 

26.8 

13 

9 

97.5 

89.4 

21 .91 

11* 

10 

101.1 

109.5 

20.63 

15 

10 

101 .6 

112.6 

11.84 

16 

10 

103.6 

96.8 

15.87 

17 

11 

106.6 

125.5 

I6.87 

18 

9 

106.6 

103.2 

15.21 

19 

10 

101*. 8 

108.5 

16.78 

20 

9 

101  .9 

107.5 

11.82 

Av  1-10 

90.58 

96.84 

Av1 1-20 

101.37 

104.11 

> 

< 

N1 

O 

95.98 

100.48 

Starts  at  042801  ,  S3 -4  Rev  445 ,  1550  A  photometer. 
320  samples,  4  counts  per  bin. 
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would  indicate  a  search  for  periodicities  that  might  be  explainable  as 
geophysical  effects  on  the  mixing  ratios  of  species  that  emit  (and 
absorb)  Lyman-B i rge-Hopf i el d  band  ultraviolet  excited  by  sunlight  between 
about  135  and  170  km  (the  satellite's)  altitude. 

In  the  absence  of  such  lower-thermospher ic  perturbat ions  the 
numbers  of  counts  would  be  expected  to  be  Po i s son-d i s t r i buted .  Student's 
t-test  shows  that  the  count-rate  data  sample  selected,  when  their  long¬ 
term  trend  is  removed  by  high  pass  filtering,  do  not  fit  this  distribu¬ 
tion.  The  Chi-square  test  shows  the  0.01  sec  radiance  samples  in  in¬ 
dividual  3-2  sec,  ^25  km-longitude  segments  to  be  normally  distributed, 
but  with  means  somewhat  higher  than  expected.  This  suggests  again  that 
short-period  (<  25  km)  atmosphere  processes  are  affecting  the  VUV 
column  emission.  As  noted  above,  periodicities  in  this  variation  can 
be  assessed  by  applying  autocorrelation  -Fourier  transform  methods  to 
the  data  stream.  We  judge  that  it  would  also  be  useful  to  investigate 
the  statistics  of  the  fluctuations  when  the  magnetic  disturbance  index 
Kp  is  higher  than  in  the  data  set  assessed  here  (2+)  ,  as  the  high- 
latitude  density  variations  are  associated  with  magnetic  activity  (Ref  **2). 

We  point  out  again  that  a  finding  of  fluctuations  in  the  1*485  - 
1650  A  FWHM  vert ica 1 -col umn  radiance  of  the  sunlit  atmosphere  means 
that  infrared  emission  when  the  M35  -  180  km  altitude  region  is 
uniformly  dosed  by  nuclear  radiations  will  be  by  and  large  similarly 
structured.  For  example  the  infrared  volume  emission  rates  from  NO 
(largely  from  0  atom  interchange  reactions)  and  0^  (thermal  and  chemi¬ 
luminescent  emission,  and  scattering  of  earthshine)  would  follow  the 
concentration  of  0  atoms;  and  the  CO^  concentration  is  expected  to  be 
coupled  to  that  of  the  heavier  minor  species.  (Conversely,  weak  clutter 
in  the  sky's  L-B-H  radiance  signal  would  mean  high  s ignal /background 
noise  in  surveillance  systems  operating  in  this  (broad)  band  of  VUV 
wavelengths.) 
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J.  Workman 

ATTN 

A.  Aikin 

Boston  College 

NASA 

ATTN 

Dept  of  Chemistry,  D.  McFadden 

ATTN 

W.  Roberts 

ATTN 

Chairman  Dept  of  Physics 

ATTN 

Science  Library,  F.  McElroy 

NASA 

ATTN 

E.  Schmerling 

California 

Institute  of  Technology 

ATTN 

R.  Schiffer 

ATTN 

S.  Trajmar 

ATTN 

N.  Roman 

ATTN 

V.  Anicich 

NASA 

Chem  Data  Research 

ATTN 

Code  JM6,  Technical  Library 

ATTN 

K.  Schofield 

NASA 

Cornell  University 

ATTN 

G.  Poppoff 

ATTN 

M.  Kelly 

ATTN 

W.  Starr 

ATTN 

N-245-3,  R.  Whitten 

University 

of  Denver 

ATTN 

Security  Officer  for  B.  Van  Zyl 

National  Science  Foundation 

ATTN 

Sec  Officer  for  D.  Murcray 

ATTN:  Div  of  Atmos  Sci 

ATTN 

R.  Sinclair 

Epsilon  Labs,  Inc 

ATTN 

C.  Accardo 

OTHER 

- 

ESL,  Inc 

University 

of  Minnesota 

ATTN 

W.  Bell 

ATTN 

J.  Winkler 

General  Electric  Co 

ATTN:  P.  Zavitsanos 

ATTN:  Technical  Information  Center 

ATTN:  M.  Linevsky 

General  Electric  Co 

ATTN:  J.  Schroeder 

General  Research  Corp 
ATTN:  J.  Ise,  Jr 

General  Research  Corp 

ATTN:  T.  Zakrzewski 

Horizons  Technology,  Inc 
ATTN:  R.  Kruger 

HSS,  Inc 

ATTN:  M.  Shuler 
ATTN:  D.  Hansen 

Information  Science,  Inc 
ATTN:  W.  Dudziak 

Institute  for  Defense  Analyses 
ATTN:  E.  Bauer 
ATTN:  H.  Wolfhard 

Ion  Physics  Corp 

ATTN:  C.  Hauer 

Johns  Hopkins  University 
ATTN:  J.  Kaufman 

Kaman  Sciences  Corp 
ATTN:  W.  Rich 
ATTN:  D.  Foxwell 

Kaman  Sciences  Corp 

ATTN:  T.  Stephens 

Kaman  Tempo 

ATTN:  DAS  I AC 
ATTN:  B.  Gambill 
ATTN:  J.  Thompson 
ATTN:  W.  Knapp 

KMS  Fusion,  Inc 

ATTN:  Library 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  D.  Divis 

Lockheed  Missiles  &  Space  Co,  Inc 
ATTN:  R.  Gunton 
ATTN:  T.  James 
ATTN:  8.  McCormac 
ATTN:  M.  Walt 
ATTN:  J.  Reagan 
ATTN:  J.  Evans 
ATTN:  J.  Kumer 
ATTN:  R.  Sears 

University  of  Lowell 
ATTN:  G.  Best 

University  of  Maryland 

ATTN:  Chemistry  Dept,  J.  Vanderslice 


University  of  Massachusetts 
ATTN:  H.  Sakai 

University  of  Minnesota 
ATTN:  M.  Hirsch 

Mission  Research  Corp 
ATTN:  D.  Archer 
ATTN:  R.  Hendrick 
ATTN:  M.  Scheibe 
ATTN:  R.  Kilb 
ATTN:  D.  Sappenfield 
ATTN:  M.  Messier 
ATTN:  W.  White 

National  Academy  of  Sciences 
ATTN:  J.  Sievers 

ATTN:  National  Materials  Advisory  Board 

Nichols  Research  Corp,  Inc 
ATTN:  N.  Byrn 

Pacific-Sierra  Research  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 
ATTN:  E.  Field,  Jr 

Pennsylvania  State  University 
ATTN:  L.  Hale 
ATTN:  J.  Nisbet 

Photometries,  Inc 
4  cy  ATTN:  I.  Kofsky 

4  cy  ATTN:  D.  Villanucci 

4  cy  ATTN:  M.  Chamberlain 

4  cy  ATTN:  C.  Trowbridge 

Physical  Science  Lab 

ATTN:  W.  Berning 

Physical  Sciences,  Inc 
ATTN:  K.  Wray 

ATTN:  R.  Taylor 

ATTN:  G.  Caledonia 

University  of  the  Commonwealth,  Pittsburgh 
ATTN:  M.  Biondi 

ATTN:  F.  Kaufman 

ATTN:  W.  Fite 

Quantum  Systems,  Inc 

ATTN:  S.  Ormonde 

R  &  D  Associates 

ATTN:  F.  Gilmore 

ATTN:  B.  Gabbard 

ATTN:  R.  Lelevier 

ATTN:  R.  Turco 

ATTN:  P.  Haas 

ATTN:  H.  Ory 

R  &  D  Associates 

ATTN:  B.  Yoon 

ATTN:  J.  Rosengren 

Science  Applications,  Inc 
ATTN:  D.  Hamlin 

Science  Applications,  Inc 
ATTN:  N.  Byrn 
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Radiation  Research  Associates,  Inc 
ATTN:  N.  Schaeffer 

Science  Applications,  Inc 
ATTN:  R.  Johnston 

Sechiist,  Professor  Chalmers  F. 
ATTN:  S.  Bowhill 
ATTN:  C.  Sechrist 

Spectral  Sciences,  Inc 
ATTN:  F.  Bein 

SRI  International 

ATTN:  R.  leadabrand 
ATTN:  R.  Hake,  Jr. 

SRI  International 

ATTN:  C.  Hulbert 

Sylvania  Systems  Group 

ATTN:  J.  Concordia 

Strategic  Systems  Div 

ATTN:  I.  Kohl  berg 

Technology  International  Corp 
ATTN:  W.  Boquist 

University  of  Texas  System 
ATTN:  J.  Browne 

TRW  Elec  ronics  &  Defense  Sector 
A i i N:  Technical  Information 
ATTN:  J.  Frichtenicht 


Utah  State  University 

ATTN:  Sec  Control  Ofc  for  K.  Baker, 

Dir  Atmos  &  Space  Sci 
ATTN:  Sec  Control  Ofc  for  D.  Burt 
ATTN:  Sec  Control  Ofc  for  A.  Steed 

University  of  Virginia 
ATTN:  R.  Ritter 
ATTN:  H.  Kelly 
ATTN:  R.  McKnight 

Vi  si  dyne,  Inc 

ATTN:  C.  Humphrey 
ATTN:  T.  Degges 
ATTN:  H.  Smith 
ATTN:  J.  Carpenter 
ATTN:  0.  Manley 

Wayne  State  University 
ATTN:  P.  Rol 
ATTN:  R.  Kummler 

Westinghouse  Electric  Corp 
ATTN:  P.  Chantry 

William  Marsh  Rice  University 
ATTN:  R.  Stebbings 

William  Marsh  Rice  University 

ATTN:  Industrial  Sec  Super,  J.  Chamberlain 
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